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SUMMARY 


This  Final  Technical  Report  presents  the  results  of  a  feasibility  investigation  of  a 
technique  for  calculating  the  far-field  radar  cross  section  of  an  object  based  upon 
measurements  made  in  the  near-field  of  the  object.  This  technique  is  an  extension  of 
existing  planar  near-field  antenna  measurement  technology,  and  is  capable  of  measuring 
monostatic  radar  cross  section  and  bistatic  radar  cross  section  at  both  narrow  and  wide 
angles. 

Included  are  the  detailed  formulation  of  the  theory  of  near-field  planar  bistatic 
radar  cross  section  measurement,  and  discussion  of  details  of  the  mechanical  scanner  and 
software  implementation.  The  comparison  of  measurement  with  predictions  is  presented; 
the  agreement  is  excellent,  and  suggests  that  a  larger-scale  demonstration  would  be 
appropriate.  Also  included  are  concepts  for  reducing  the  amount  of  data  required  for 
reconstructing  radar  cross  section,  and  a  discussion  of  limitations  of  this  method  of  radar 
cross-section  measurement. 

Specific  recommendations  are  presented  for  technology  development  areas  that 
should  be  pursued  to  mature  this  measurement  technique  into  a  viable,  operational 
technology.  Among  those  areas  are  calibration,  data  handling,  computational 


FOREWORD 

The  BDM  Corporation,  1801  Randolph  Road  SE,  Albuquerque,  NM  87106,  is 
pleased  to  submit  this  report,  titled  "Final  Technical  Report  for  Near-Field  Bistatic  RCS 
Measurement,"  to  the  Rome  Air  Development  Center  as  required  by  CDRL 
DI-A-3591A/M. 

This  document  presents  a  description  of  the  work  performed  under  contract  number 
F19628-86-C-0208  during  the  period  of  September  30, 1986  to  March  20, 1989. 
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CHAPTER  I 
INTRODUCTION 

Radar  cross-section  (RCS)  measurement  is  an  important  ingredient  of  design 
verification  and  maintenance  of  contemporary  aircraft  and  missile  systems.  RCS  is 
becoming  more  and  more  important  for  both  tactical  and  strategic  weapons  systems,  due  in 
part  to  increasingly  sophisticated  radar  systems  and  the  resulting  need  for  low-observable 
aircraft  and  missiles. 

BDM  is  currently  under  contract  to  Rome  Air  Development  Center  to  validate  a 
technique  for  calculating  far-field  monostatic  and  bistatic  RCS  using  planar  measurements 
in  the  near-field  of  a  test  object.  This  technique,  developed  by  BDM  in  1984,  is  of 
increasing  relevance  to  the  national  need  in  this  area,  particularly  in  view  of  some  of  the 
unique  capabilities  that  the  planar  near-field  technique  may  provide. 

The  goal  of  thi3  contract  effort  is  to  demonstrate  the  feasibility  of  near-field 
measurement  of  bistatic  radar  cross  section.  The  proposed  technical  approach  is  divided 
into  six  tasks  which  may  be  summarized  as  Mows: 

.  Task  1:  Formulation  of  probe-corrected  near-field  RCS  theory 

Task  2:  Development  of  software  to  compute  scattering 

Task  3:  Design  and  construction  of  scanner 

Task  4:  Comparison  of  measurements  with  solutions 

Task  5:  Data  reduction  investigation 

Task  6:  Investigation  of  method's  limitations 

Those  six  tasks  are  addressed  in  Chapters  II  through  VII  of  this  final  report,  while 
Conclusions  are  presented  in  Chapter  VIII. 

It  is  important  to  realize  that  some  of  today's  RCS  measurement  requirements 
simply  cannot  be  met  with  existing  measurement  technologies.  For  example,  measuring 
the  RCS  of  low-observable  aircraft  is  difficult  on  conventional  far-field  ranges:  the 
signature  is  small  and  hard  to  measure;  range  effects  (ground  and  air  scattering)  degrade 
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accuracy  and  repeatability;  airborne  surveillance  during  long  observation  times 
compromises  security;  and  bistatic  RCSis  difficult  to  measure  (particularly  at  small 
angles). 

BDM  has  been  working  on  near-field  RCS  measurement  and  prediction  concepts 
since  1983.  la  1984  Mike  Dinallo  developed  and  published  the  mathematical  foundations  of 
a  proposed  neat-field  RCS  (NFRCS)  measurement  approach  which  we  called  "planar 
near-field  RCS  measurement."  In  1987  Rogers  and  Farr  published  the  explicit  solution  to 
Dinallo's  scattering  equations  and  described  BDM's  ongoing  work  in  near-field  RCS 
measurement. 

Our  planar  NFRCS  measurement  technology  is  an  outgrowth  and  extension  of  the 
near-field  antenna  measurement  theory  and  techniques  developed  by  the  National  Bureau 
of  Standards  (NBS)  in  the  1960s.  The  N3S  theoretical  anti  experimental  programs 
demonstrated  that  traditional  far-field  antenna  patterns  can  be  calculated  based  upon 
pattern  data  measured  in  the  near-field  of  an  antenna.  These  results  were  backed  up  by 
extensive  error  analyzes  and  validation  tests  of  antennas  on  near-field  and  far-field  ranges. 
The  analyses  and  tests  showed  that  patterns  measured  using  near-field  techniques  are 
more  accurate,  repeatable,  and  economical  to  obtain  than  patterns  measured  on  traditional 
antenna  ranges.  The  near-field  measurement  facilities  are  also  smaller  than  far-field 
ranges  and  are  fully  enclosed,  providing  excellent  security  for  sensitive  military 
applications.  The  NBS  near-field  antenna  measurement  methods  are  now  in  daily  use  by 
most  major  antenna  fabrication  and  test  facilities. 

Since  1984  BDM  has  been  extending  the  NBS  near-field  antenna  measurement 
theory  to  include  the  near-field  RCS  measurement  problem.  Our  work  has  been  both 
theoretical  and  experimental;  we  have  developed  the  mathematics  of  near-field  RCS 
measurement  and  have  carried  on  an  active  experimental  program  to  validate  the 
mathematics.  Our  planar  near-field  RCS  measurement  technique  is  the  result  of  this 
work. 
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The  key  feature  of  this  technique  is  the  mathematical  algorithm  that  allows  us  to 
efficiently  compute  the  RCS  of  an  object  based  upon  many  scattering  measurements  made 
near  the  object.  The  measurements  are  conceptually  easy  to  obtain  by  a  computerized 
data  control  and  acquisition  system  similar  to  the  one  we  have  prototyped  in  the  BDM 
Laboratory. 

The  unknown  object  is  illuminated  by  a  single  broad-beam  transmitting  antenna 
and  the  scattered  signal  from  the  object  is  received  by  one  or  more  receiving  antennas.  The 
transmitting  and  receiving  antennas  are  in  the  near-field  of  the  object  (typically  within  a 
few  feet  of  it  when  using  gigahertz  frequencies).  A  computerized  data  acquisition  system 
controls  where  the  transmit  and  receiving  antennas  are  placed  (using  servo  control)  and 
measures  and  stores  the  received  signals.  The  transmitting  and  receiving  antennas  are 
moved  around  in  a  plane  (i.e.,  a  planar  scan  pattern),  so  our  measurement  approach  is 
more  correctly  called  near-field  bistatic  RCS  measurement  using  planar  scanning. 

Other  noteworthy  features  are  that  bistatic  measurements  are  feasible  both  at  large 
and  small  angles,  monostatic  measurements  are  feasible,  sensitivity  is  excellent  (so 
low-observable  can  be  measured),  and  the  measurement  facility  is  totally  enclosed  (which 
enhances  security). 
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CHAPTER  H 

FORMULATION  OF  NEAR-FIELD  THEORY 
A.  INTRODUCTION 

In  this  chapter  we  present  the  analysis  of  near— field  bistatic  scattering  data.  This 
analysis  is  based  upon  Dinallo's  [16]  formulation  of  bistatic  scattering  in  terms  of  the  plane 
wave  scattering  matrix.  Note  that  for  the  built  of  this  section  we  carry  out  the 
mathematics  in  the  unprimed  coordinates  of  the  object  under  test. 
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SCATTERING  MATRIX 

The  goal  is  to  calculate  fcr  an  arbitrary  object  a  scattering  matrix 

hitfih  w^>. 

which  describes  the  scattering  of  an  incident  plane  wave  by  that  object,  and 
h  =  kxex  +  kyey  +  kz§z 
1  =  lx«x  +  lySy  +  lzez. 

It  is  worth  noting  explicitly  that  the  dements  of  the  scattering  matrix  l11(k,l)  are 
specified  by  means  of  two  polarization-related  indices  and  two  propagation  vectors.  For 
example,  Ijj^k,!)  refers  to  the  ^-component  of  the  angular  spectrum  of  the  wave  that  is 
scattered  when  the  object  is  illuminated  by  the  ^-component  of  the  angular  spectrum  of 
the  incident  wave.  The  incident  wave  has  propagation  vector  T,  and  the  scattered  wave 
has  propagation  vector  k.  Taking  all  four  elements  together,  Ij^k,!)  specifies  how  the  8 
and  #  components  of  the  angular  spectrum  of  the  incident  wave  1  are  scattered  into  8  and 
<t>  components  of  the  angular  spectrum  of  a  scattered  wave  k. 

The  incident  and  scattered  waves  are  generally  a  superposition  of  plane  waves  which 
we  will  modd  as  a  continuum  of  plane  waves.  A  complicated  wavefront  illuminating  the 
object  may  be  decomposed  into  a  spectrum  of  plane  waves.  Measuring  the  transverse 
components  E^  and  E.y  of  that  E-fidd  in  some  plane,  such  as  z  =  zQ>  the  angular 


n*i 


spectrum  (of  plane  waves)  is  the  Fourier  transform  in  the  x-y  plane 
- 1 1  E,(r')  e  dr', 


r'=r  e  +re  +  z  e  . 
xx  77  o  z 

The  multitude  of  plane  waves  that  makes  up  that  angular  spectrum  of  Incident  plane  waves 
is  scattered  by  the  object,  yielding  an  infinity  of  emergent  plane  waves  whose  angular 
spectrum  is 


i2(5)= yhhg+ty&ht 

[WM>  *u**&^]  [ 


^,1)1  IU(1)  • 

Note  that  the  0-  and  ^-components  of  the  angular  spectrum  of  the  scattered  plane 
waves  are  a  linear  combination  of  both  the  0-  and  the  ^-components  of  the  angular 
spectrum  of  the  incident  plane  waves. 

For  a  source  at  having  angular  spectrum  1^(1),  the  scattered  angular  spectrum  at 
the  origin  is  denoted  by  F(rltk),  where 


F (?jA- 


F/rltS) 

*11 0^’^ 


*10^  -l!-f, 

6 


j  j  I  \  [  h/t)  j 

(As  written  here,  l'(r1Jc)  is  related  directly  to  the  sampling  plane  coordinates  rather  than 
the  object  coordinates.)  The  f&r-field  scattered  £  at  ?2  is  calculated  with  the  results  from 
a  previous  section.  Hence, 


.  -2rikcosg 


lHU) 


where  k  is  parallel  to  ?2.  If  the  incident  wave  is  in  fact  planar,  then  the  angular  spectrum 


of  the  incident  wave  is 

Pi*] 

IjO)  = 

IV 

and  the  integral  above  reduces  to  a  matrix  multiplication.  The  far-field  scattered  E-field 
is 


Given  illumination  of  the  test  object  by  the  transmit  antenna  in  transmit 
orientation  #1,  we  will  first  calculate  the  9-  and  ^-components  of  the  scattered  angular 
spectrum  over  the  area  scanned  by  the  receive  probe  antenna  in  order  to  calculate  the 
scattering  matrix.  Next,  the  transmit  probe  antenna  is  rotated  90°  to  illuminate  the  test 
object  with  a  different  angular  spectrum  of  plane  waves,  and  we  then  repeat  the 
measurements  and  calculations  in  order  to  determine  the  6-  and  ^-components  of  the 
scattered  angular  spectrum  for  this  second  transmit  probe  antenna  orientation.  From  these 
two  sets  of  data,  we  can  calculate  the  manner  in  which  the  test  object  scatters  the  &-  and 
^-components  of  the  incident  wave  into  the  9~  and  ^-components  of  scattered  waves. 

The  laboratory  measurements  for  bistatic  near-field  scattering  consist  of  gain  and 
phase  measurements  made  by  a  receive  probe  as  it  is  swept  through  a  pattern  of  probe 
locations.  The  receive  probe  scan  pattern  is  swept  repeatedly  (and  gain  and  phase 
measurements  made)  for  a  set  of  transmit  probe  locations,  as  the  transmit  probe  itself  steps 
through  a  set  of  locations  in  the  transmit  probe  scan  pattern. 

The  transmit  probe  and  receive  probe  we  use  are  identical,  although  they  need  not 
be  so.  The  receive  probe's  receiving  characteristic 


and  the  transmit  probe's  transmitting  characteristic 

rw1) 

10^)=  i  (!)  ’ 

are  assumed  to  be  known  (either  analytically  or  empirically).  As  in  the  case  of  near-field 
antenna  measurements,  neither  probe  antenna  should  have  nulls  in  a  direction  r  for  which 
one  wishes  to  calculate  the  scattering  matrix. 


In  this  section  we  outline  the  approach  we  take  to  determine  the  scattering  matrix 

ljj(k,!)  The  scattering  matrix  l^k.l)  of  an  object  is  a  tensor  of  rank  two  that 

describes  how  an  incident  plane  wave  (1)  with  arbitrary  polarization  is  scattered  by  that 

object  and  transformed  into  one  or  more  scattered  plane  waves  (k).  The  scattering  of  an 

arbitrary  wavefront  follows  directly,  since  the  superposition  principle  allows  us  to  create 

complex  wavefronts  as  a  superposition  of  plane  waves. 

The  scattering  object  is  at  the  origin,  and  we  deSne 

!  =  propagation  vector  of  incident  wave  (from  transmit  probe), 

* 

Tj  =  coordinates  of  source  of  incident  wave  (transmit  probe), 

* 

k  =  propagation  vector  of  scattered  wave  (to  receive  probe),  and 
* 

r0  =  coordinates  of  receive  probe. 

The  measurement  equation  for  near-field  bistatic  RCS  is  [17] 

V  V2)  =  {[$!<*)•  {}  in(5J)-i10(T)  e"i!'ri  d£  e+ik'r2dKdL. 

{5}  {3}  {2}  {  1}  {4} 

The  scalar  bQ  is  a  measurable  quantity  that  is  the  complex  gain  factor  describing 

the  transmission  path  from  the  transmit  probe  to  the  receive  probe  via  the  scattering 
object.  The  incident  wave  is  represented  by  the  factor  {2}.  The  exponential  {1}  represents 
propagation  of  the  incident  wave  from  the  transmitting  antenna  to  the  object,  and  the 
multiplication  with  {3}  represents  scattering  by  the  object.  The  resulting  scattered 


spectrum  of  waves  propagates  {4}  from  the  object  to  the  receive  probe;  {5}  represents  the 
interaction  of  the  propagated,  scattered  wave  with  the  receiving  probe  itself.  The  gain  and 
phase  ratios  measured  by  the  Hewlett-Packard  network  analyzer  are  the  values  of  bQ. 

If  we  consider  the  scattering  object  as  a  transmitting  "antenna,"  the  complicated 
wavefront  generated  by  it  can  be  represented  by  a  superposition  of  plane  waves.  We 
denote  that  superposition  by  f^(r^jc),  where  r^  is  the  location  of  the  actual  transmit 
source  and  k  is  the  wavevector  of  the  scattered  plane  waves.  Since  the  scattered  spectrum 
is 

4 

*!(?!■*)-  !n(k,i)-i10(i)  'fldL, 

the  measurement  equation  above  becomes 

„  .  ff.  .......  +ik-r0 

b0(rrr2>“  I01(k)*F1(r1^)  e  2dK. 

This  is  the  measurement  equation  for  near-field  antenna  measurement  which  may  be 
solved  to  determine  P^r^k). 

Furthermore,  the  above  two  equations  are  very  similar,  except  that  b^r^,^)  is  a 
scalar  equation  while  Fj(rj,k)  is  a  vector  equation.  The  techniques  of  an  earlier  section 
can  be  applied  to  solving  for  ^(r^jc),  except  that  two  independent  solutions  must  be 
found  (representing  the  eg  and  e^  components  of  ^(r^k)). 


Using  the  above  definition  of  the  scattered  spectrum  f(?j,k),  the  measurement 
equation  can  be  written 

boi(?i*?2> =  ioiOO-W) e+  *  2(& 
which  is 'a  Fourier  integral  that  can  be  inverted  formally  [18]  as 


The  integral  is  over  ail  positions  of  the  receive  probe,  holding  the  transmit  probe 
position  fixed.  Since  the  receive  probe  moves  in  the  plane  z  —  d  (a  constant),  we  define 

®2  *  r2xex  +  r2y®y  +  ^‘®z’ 
and  the  integral  becomes 

-idk  rr  .A  * 

«  ^p-  b0l(^,?2) e 

We  define  the  right-hand  side  of  the  equation  as  a  coupling  product  D,  (r^lc),  so 

10;(5)-W)  =  Dl^lA- 

where  D,  (?,,£)  is  the  calculable  quantity 

-  -  ~idkz  rr  „  ,  +SK-L  . 

Dl(rl,k)  *  ” 2  bol^rl,r2^ e 

=wM/ri*)+w{S)W'i*)- 

Let  the  measurements  be  repeated  with  the  same  transmit  probe  but  with  a 
different  receive  probe  antenna  that  has  receiving  characteristic 
161*00 

.W£)[ 

obtained  by  rotating  the  receive  probe  by  90?  Calling  the  new  set  of  measurements 

bo2^rl’r2^’  we  can  compute  another  coupling  product  D2(rj,k)  as 

,  ,  ~idiz  ff  „  ,  -iK-jL  , 

°2(rl’k)  =  -pi  1  ibo2<rl’r2^  e  ^  d&2’ 

*  I01^)F1^1>5)  + 

If  the  determinant 

i3  non-sero,  then 


*01® 


n-6 


1 

In  summary,  the  transmit  probe  antenna  is  in  an  orientation  that  we  define  as  "TX 
orientation  #1.H  For  each  transmit  probe  location  r^,  we  measure  an  array  of  scalar  values 
as  the  receive  probe  (in  RX  orientation  #1}  is  moved  through  its  scan  pattern.  A  second 
array  of  scalar  values  is  measured  with  the  receive  probe  rotated  by  90°  (RX  orientation 
#2).  These  two  arrays  of  measured  data,  called  b^r^)  and  bo2(r1,r2),  are  used  to 
calculate  Fj(i'pk).  Fj(r^,k)  is  a  vector  quantity  that  is  defined  for  every  transmit  probe 
location  and  every  possible  k. 

F-  gQMPUTINgialMl 

In  order  to  calculate  the  scattering  matrix,  it  is  necessary  to  determine  a  second, 
independent  scattered  spectrum  ?2(rj,k)  that  represents  Use  scattered  spectrum  with  the 
object  illuminated  by  a  different  source  spectrum.  Since  the  transmit  probe  is  linearly 
polarized,  we  can  create  a  different  source  spectrum  by  rotating  the  transmit  probe 
antenna. 


With  the  transmit  probe  antenna  rotated  by  90°  (TX  orientation  #2),  we  repeat 
the  entire  series  of  measurements  and  calculations  described  in  the  previous  section. 
For  each  transmit  probe  location  r,(  we  measure  an  array  of  scalar  values  as  the  receive 
probe  (in  RX  orientation  #1)  i3  moved  through  its  scan  pattern,  A  second  array  of  scalar 
values  is  measured  with  the  receive  probe  in  RX  orientation  #2.  These  two  arrays  of 
measured  data,  called  bo3(r1,r2)  and  bo4(rl,r2),  are  used  to  calculate  F^.k). 

The  pertinent  equations  are  summarized  below: 


b°3(?l’?2)  -  J 


I^-F/r^e""  ‘2dk, 


+ik-r„ 


i(k)*f2(ri$  *  pr  bo3^i>^e  *  2<ih' 
eidk*  ff  ,  .  -dfc.fi,  A  “ 


2  Jtyvk) 


or,  explicitly, 


!oiCHM  -  hnl&fy/1!*) + 


idk„  pj-  .  *  ^ 

«  ~jr —  bo4(r1,r2)  e  '  dSj 
=  I6l^^F2^rl’^)  •*■  r6i^^F2^?iA- 


^2'k) 2  -I6i^k)Ioi^) 

is  non-zero,  then 


Pj^Jc)  -  -D3(ri^r6i^;x  ’ 

AjCr^k) 


F,  ,(r, ,k)  =  .  ~ds(ti^)Iqx^?i^) 

2* 1  ~ 


Refemag  to  the  general  bistable  Jneasnieosat  equation 

b0(rl»r2)*  foi^)-[{{ in&MJ0(I) «  ^  ridl]e+'JC  \M, 

we  have  calculated  f^i)  and  fyxji), 

-  If  iu&Htf!)  e^1  d£ 


f/ri’k)=|  iu(k.I)-Jio(t)e‘,i'rid£. 
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are  Fourier  integrals  that  my  be  inverted  to  give 
i  /t\ _ I _ fCrrl!  fl 


An(k,I).!10(I)  =  Apj  j  F^.I)  e  fl  dtj 

““  ~J~  ^2^rl’^e  ^  dr  j. 


This  time  tbe  integrals  are  over  all  positions  of  the  transmit  probe,  holding  ibe 
receive  probe  position  fixed.  Since  the  transmit  probe  moves  in  the  plane  z  =•  d,  iet 


Rl  =  rlx®x  +  Vy  +  d'V 
and  the  integrals  become 

+idl,,  fr  +j£. 

“Mdl,.  rr 


+i£.H1 


in(k,l)4io(l)»S_r..|lf2(r1.l)e  ‘dRj. 

We  define  a  pair  of  vector  coupling  products  ($(kj)  and  $'(»,!)  as 


so  that  we  can  write 


ln^)-Uh  = 

ij,(h,I)*?jp{T)  *  3'(k,T). 
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Expanding  the  above  equation  in  components  and  dropping  the  explicit  (k,I) 


dependence, 


W  Ji0^  0  0  |  I1!!**!  j'  Q$ 


0  0  hw  ho#  .  hl0ti 


Hit  Hot  0  0 


0  0  hoffhotj  [^wj  [  %  j 
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wLiifi  has  solutions: 


t  ~  Qfto* 

llW  -  A  — 


hiu  - - 2S"-~- 

provided  the  determinant 

A3l10^10^-I105I10^ 

is  ncaaao  for  particular  (k,l).  This  completes  the  calculation  of  the  scattering  matrix 

hlW)- 

Note  that  all  of  the  above  calculations  are  in  the  coordinate  system  of  the  object 
under  test.  When  making  the  actual  measurements,  one  must  immediately  convert  from 
the  coordinates  system  of  the  measurement  table  to  that  of  the  object  under  test. 


IMO 


CHAPTER  III 

DEVELOPMENT  OF  SOFTWARE  TO  COMPUTE  SCATTERING 


A.  INTRODUCTION 


The  software  to  perform  the  scattering  calculations  is  a  straightforward 
implementation  of  the  algorithm  presented  in  the  previous  chapter.  The  heart  of  the 
software  is  the  two-dimensional  fast— Fourier  transform  subroutine;  a  number  of  such 
subroutines  are  available  in  the  open  literature. 


There  are  a  number  of  software  issues  stemming  from  the  discretization  that  is 
implicitly  employed  to  calculate  discrete  equivalents  to  continuous  integrals.  These  issues 
are  discussed  iu  the  following  sections. 

B-  BET  ORIGIN  AND.  PHASE  SRIFT 

In  this  section  we  discuss  sampled-data  evaluation  of  Fourier  integrals  by  the  FFT 
algorithm.  For  rigorous  discussions  refer  to  the  literature  [12J. 

Numerical  evaluation  of  integrals  of  the  form 


H  («)< 


h(t)  elut  dt 


is  feasible  for  bandiimited  h(t)  using  sampled  data  and  discrete  Foi’rier  transform 
techniques.  The  continuous  function  h(t)  is  sampled  at  intervals  of  St  chosen  by  Nyquist's 
criterion.  The  sampled  h(t)  is  a  finite  sequence 

{h'(n)  =  h(t)  j  1 3  +  (n— n»l,..  .N} 

3h'(l),h'(2),...h'(N) 

which  is  implicitly  periodic  in  n  with  period  N.  The  DFT  of  h'(n)  is 
N  — 2ri(n— llfk-ll 

K'(k)=  l  h'(n)e  S  ,  k  *  I,  ...N 
n=l- 

where  (for  even  N) 


m-i 


u  = 


(k-1 )  6oj  k  =  l ,  2, .. . — rf'-f  1 

[k— (N+l)]&n  k  =-Tr-f2,-rj-+3,  ...  N 

H(w)  =  fc-H'(k). 

For  example,  with  N=8, 

H'(l)  =  H(+0  Su) 

H'(2) »  H(+l5w) 

H'(3)  *  H(+2dwj 
H'(4)  =  H(+3  Su) 

H'(5)  =  H(+4M  =  H(-4&>) 

H'(6)  =  H(-35u) 

H'  7)  a  H  -2 Su) 

H'(8)  =  H(-lH 

In  a  sampled— data  implementation,  Fourier  integrals  are  conveniently  calculated 
using  the  discrete  Fourier  transform  (DFT)  in  the  form  of  the  Fast  Fourier  Transform 
(FFT)  algorithm.  However,  there  is  a  phase  shift  implicit  in  the  DFT  that  arises  because 


the  origin  of  the  transform  is  at  the  first  point  processed. 
We  compute  the  above  integral  by  a  summation: 


where  I  ,  I  are  indices  of  a  measurement  data  array,  and  K  ,  K  are  (in  this  instance 
only)  indices  of  the  transf  <nred  data  array.  Ix  and  Kx  range  from  I  to  Nx,  and  Iy  and  Ky 
range  from  1  to  Ny. 

Idealizing  the  sampling  process  somewhat,  the  sampled  function  f  (L,I  )  at  given 

“  *  J 

(Is,Iy)  corresponds  to  the  continuous  function  f(xj)  at  the  point 

yyy)  *  ^xo+Px-^x-  v!Irl3V- 

nr-2 


so 


The  measurement  data  array  f3(-)  is  defined  over  some  set  of  (Ix,Iy)  that  does  not 
necessarily  include  the  coordinate  system  origin.  In  order  that  the  computed  phase  be 
referenced  to  that  origin,  a  phase  factor  must  be  included  to  specify  the  phase  of  the  (1,1) 
point  relative  to  that  origin: 


Nx  Ny 

I  I  fM),Kp 


lx=Uy=l 
where 


±2ri- 


•exp 


'  (K  -1 )  (I  -1) 
7  N;;y— 


#Kx,Ky)  =  exp 


*2*i- 


.(K -l)o 


•exp 


’  (Kv-l)a 

*2ri . -V-X 


and 


ax~rox^x, 

ay  ~  roy^V 

No  modification  to  the  FFT  algorithm  is  required,  since  the  phase  factor  ^Kx,Ky) 
can  be  calculated  and  multiplied  into  the  transformed  data  array  generated  by  the  FFT 
algorithm. 

C.  DISCUSSION  OF  CALCULATED  SCATTEEING.MAIRK 

For  a  given  test  object,  our  laboratory  measurements  and  data  analysis  produce  a 
set  of  scattering  matrices  of  the  form 

‘llxx^i’V  ^lxy^i’ty 

WVl>  WVj> 

where  each  matrix  consists  of  four  complex  scalars  (Illxx,  Ijlyx,  ^jyy) tliat  define 
the  scattering  matrix  for  given  values  of  K  and  L.  As  before,  L  is  the  projection  of  the 
incident  (illuminating)  plane  wave  propagation  vector  1  onto  the  x-y  plane,  and  K  i3  the 
projection  of  the  scattered  plane  wave  propagation  vector  k  onto  the  x-y  plane. 

Since  the  discrete  Fourier  transform  (DFT)  is  used  in  the  sampled-data 
implementation  of  the  integrals,  K  and  L  can  take  on  only  a  discrete  set  of  values  that  are 


ra-3 


determined  by  the  measurement  grid  size,  the  spatial  sampling  interval  &,  and  the  "rules" 
of  the  DFT.  Fourier  interpolation  can  be  used  to  increase  the  resolution  of  the  grid  of 
values  of  K  and  L 


Suppose  the  transmit  probe  scan  pattern  is  a  square  consisting  of  points  on  a 
side,  with  spatial  sample  interval  &.  The  propagation  vector  of  the  transmitted  wave  (i.e., 
the  plane  wave  that  is  incident  upon  the  test  object)  is  denoted  by 

UlA  +  1y«y  +  1A* 
and  its  projection  onto  the  x-y  plane  is 

Lsi  e  + 1  e  . 
xx  y  y 

Then  the  x—  and  the  y-componenis  of  I  take  on  the  discrete  values 

f.  NtJ 

lx  =  i - 5-  i  =  0,2, ...N^ 

*wJ-°A~Ntx 

where 

d.  _  2  x 

*■“  TS* 

is  the  spacing  between  adjacent  samples  in  "k-space".  The  maximum  magnitude  that  »x 
or  ly  can  have  is 

max  =  ST' 

However,  and  ly  are  constrained  by  the  additional  requirement  that  lz  be  real, 


x  y  z  L  A  J’ 

and  propagating  waves  correspond  to  those  1  for  which  lg  is  real.  If  lx  or  ly  gets  too  large, 
lz  becomes  imaginary,  so  we  require 
i2  .  A  ,  f  2x  I2 

‘x  +  1y-  iTi* 


Since  A-1  is  the  highest  spatial  frequency  that  can  be  present  due  to  propagating 
waves,  it  is  the  Nyquist  frequency  [13]  of  the  signal  to  be  sampled,  and  the  required 
sampling  interval  is  A/2.  If  the  sample  interval  is  indeed  S%  =  A/2,  the  domain  of  valid 
(1^  J  )  is  a  circle  in  the  (i  j)  plane  that  is  exactly  inscribed  in  the  square  defined  by 

'  J  f,  Ntx  Ntxl  .  Ntx  Ntxll 

■I*  — r-T  •»“  — T' T’jJ' 

so  a  fraction  ~  21%  of  the  computed  U  is  not  useful  to  us. 

x  xj  yj 

A  similar  situation  exists  with  regard  to  the  range  of  values  of  the  propagation 

vector  of  the  scattered  plane  wave.  Suppose  the  receive  probe  pattern  is  a  square  with 

points  on  a  side  and  spatial  sample  interval  6c.  The  propagation  vector  of  the  scattered 

wave  is  denoted  by 

k  =  k  e  +  k  e  +  k  e  , 
xcx  T  77  22 

and  its  projection  onto  the  x-y  plane  is 
^  *  *x*x  +  ky V 

The  x-  and  the  y-components  of  k  take  on  the  discrete  values 
kxj=  ‘ — T  ^rx’ ' =  °,2,-Nrx 


kXj  =  p — T  *rx’j  =  0’2’-Nr 


:"'TP“d 


■*RH‘ 


Conservative  signal  processing  technique  calls  for  using  a  sampling  interval  smaller 
than  the  A/2  dictated  by  the  Nyquist  limit.  In  the  present  instance,  the  physics  of  the 
wave  propagation  very  effectively  bandlimits  the  signal  by  imposing  a  very  sharp  cutoff  for 


spatial  frequencies  beyond  1/A.  This  i3  in  contrast  to  typical  signal  acquisition  scenarios  in 
which  a  signal  generally  has  components  above  ^Nyquist* ttere  0n3  mu3t  use  a  sharP  cut0^ 
filter  and  in  addition  sample  at  a  somewhat  higher  rate  than  the  Nyquist  theorem  requires. 
We  conclude  that  one  should  probably  use  &  <  A/2  by  perhaps  5%  or  so. 

Given  the  good  A/D  resolution  (12  bits),  adequate  floating-point  precision  and 
dynamic  range  (I.E.E.E.  standard  floating-point  format),  and  a  relatively  quiet 
measurement  location,  we  ignore  some  of  the  common  reasons  for  sampling  at  above  the 
Nyquist  frequency,  namely  quantization  and  numeric  dynamic  range. 

Since  the  experimentally-determined  scattering  matrices  are  defined  at  discrete 
values  of  K  and  L,  it  is  convenient  to  construct  one's  theoretical  models  of  scattering  such 
that  the" precise  values  of  K  and  L  calculated  by  the  analysis  software  can  be  automatically 
plugged  in  to  yield  the  theoretical  scattering  values.  From  an  algorithmic  standpoint  thi3 
corresponds  to  implementing  the  theoretical  or  numerical  model  of  scattering  as  a 
subroutine  that  has  as  input  the  values  of  K  and  L  for  which  a  theoretical  prediction  is 
needed. 

F.  OYER-SAMPLING 

Under  ideal  conditions  the  gain  and  phase  signals  from  the  network  analyzer  need  to 
be  sampled  only  at  A/2  intervals  (or  slightly  more  often  if  one  is  near  ihe  reactive 
near-field  of  the  test  object).  In  this  application  we  oversample  by  a  factor  of  ten. 
Reasons  for  this  are:  (1)  noise  generated  in  the  electronics  for  gain  and  phase  detection 
smears  out  the  spectral  content  of  the  signal  being  measured;  (2)  a  general  rule  of  thumb  in 
digitizing  and  processing  noi3y  signals  is  that  one  should  digitize  at  five  to  ten  times  the 
Nyquist  rate;  (3)  filtering  techniques  can  be  used  to  reduce  the  (uncorTelated)  noise  on  the 
signal;  and  (4)  non-linear  filtering  techniques  can  be  used  to  detect  and  correct  invalid 
phase  measurements. 

The  invalid  phase  measurements  occur  because  the  phase  detection  circuit  in  the 
network  analyzer  updates  its  output  asynchronously  with  respect  to  the  A/D  converter,  so 


it  is  possible  for  the  A/D  to  sample  the  phase  when  the  circuit  is  "wrapping"  around  from 
+180  to  -180  degrees  (or  vice  versa).  Measuring  extra  points  allows  us  to  detect  and 
correct  the  invalid  phase  values.  Spatial  filtering  to  improve  signal-to-noise  ratio  (SNR) 
is  also  practical  when  oversampling  is  performed. 

The  figures  in  this  section  were  constructed  using  measured  data  from  two  x-axis 
scans  (made  on  6/4/87)  of  10  Ghz  bistatic  RCS  from  the  6-inch  aluminum  sphere.  The 
TX  probe  was  held  stationary,  and  the  RX  probe  was  scanned  in  the  +x  direction  using  a 
200  hz  digitizing  dock  and  the  usual  probe  velocity  of  29.8  cm/sec.  Sample  interval  is 
calculated  as  0.149  cm,  corresponding  to  oversampling  by  a  factor  of  ten  at  10  Ghz. 

The  probes  were  open-end  X-band  waveguide  with  absorber  collars  to  limit 
low-angle  radiation  and  an  absorber  barrier  between  the  RX  and  TX  probes  to  limit  direct 
pro be-to-p robe  coupling.  The  TX  probe  was  driven  with  a  20-watt  (nominal) 
traveling-wave  tube  (TWT)  amplifier,  and  a  preamplifier  was  used  on  the  RX  probe. 

A  high  SNR  scan  was  obtained  using  maximum  drive  to  the  TWT  amplifier.  A  low 
SNR  scan  was  obtained  immediately  after  the  high  SNR  scan  but  with  the  TWT  drive  level 
reduced  by  15  db.  Additional  measurement  noise  was  introduced  by  the  lowered  reference 
channel  signal  at  the  network  analyzer. 

G.  ClfimZATIOKRATE 

Since  the  quantity  of  data  required  by  the  near-field  technique  is  already 
formidable,  it  is  preferable  to  store  the  minimum  number  of  values  necessary  for  the 
reconstruction  of  the  scattering  matrix.  The  signal  bang  measured  contains  no 
components  above  spatial  frequency 

^signalpnax  “  *T~  hydes/meter, 

so  digitizing  the  signal  at  a  sampling  frequency  of  2/A  samples/meter  would  theoretically 
capture  all  of  the  spectral  content  of  the  signal  For  reasons  mentioned  above,  digitizing  at 

ten  times  this  rate  is  more  appropriate,  so  we  choose: 

20 

fs  =  — j~  samples/meter. 
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Then 


&  =  fy  =  -4 —  meter/sample 
*s 
A 

=  tit 

and 

fNyquist  =  -m-cycle3/meter 

is  the  highest  spatial  frequency  that  can  be  detected  without  aliasing. 

The  bandwidth  of  the  network  analyser's  gain  and  phase  detection  circuits  is  10 
kHz,  so  samples  should  be  taken  at  intervals 

&>>irraE 

»  0.1  ms 

to  ensure  that  the  noise  in  the  samples  is  uncorrelated.  Allowing  an  order  of  magnitude 
leeway,  we  set  a  lower  bound  of  one  millisecond,  which  corresponds  to  a  spatial  interval 


&=■  v  St, 

where  v  is  the  speed  of  the  probe  antenna.  Typically  v  =  0.30  m/s,  so  &  =  .0003  meter. 

At  10  Ghz,  A/2  =  .015  meter,  so  this  &  represents  oversampling  the  data  by  a  factor  of 

A/2  _  .015  m  »„ 

~T5c - .0003  m — 50- 

which  more  than  meets  the  oversampling  goal  mentioned  above.  The  bandwidth  of  the 
receiver  is  large  enough  to  allow  much  oversampling  without  compromising  the  noise 
characteristics  of  the  sampled  data. 

Conclusions  to  be  drawn  are:  (1)  oversampling  can  be  performed  without  reducing 
the  RX  probe  scan  speed,  and  (2)  the  measured  noise  will  be  uncorrelated  from  sample  to 
sample.  Note  that  the  computer  system  is  idle  during  the  time  required  to  return  the  RX 
platform  for  the  next  scan,  so  the  Sitering  adds  no  time  penalty. 

H.  PHASE  WRAP  DETECTION 

The  phase  wrap  problem  is  difficult  to  observe  in  standard  gain  and  phase  plats 
(Figures  HI— 1  and  HI-2).  The  incorrect  measurement  occurs  when  the  phase  is  changing 
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Figure  m-l.  Gain  and  Phase  (High  SNR) 
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Figure  HI-2.  Gain  and  Phase  (Low  SJJR) 


m-10 


from  +180  to  -180  degrees  (or  vice  versa),  and  the  value  measured  is  somewhere  between 
180  and  -180,  so  that  it  appears  approximately  in  the  proper  place  on  the  plot  and  is  easily 
overlooked. 

The  invalid  sample  often  shows  up  dearly  as  a  "spike"  in  plots  of  the  real  and 
imaginary  components  (Figures  HI-3  and  IH— 4).  The  real  and  imaginary  signal 

components  are  computed  3S 
R(A,v>)  =  A  cos  ip 
I(A,y)  =  A  sin  ip 
so 

R(A,±180)  =  —  A 
I(A,±180)  =  0. 

On  plots  of  real  and  imaginary  components,  one  looks  for  phase  errors  when  both 
(1)  R  is  far  from  aero  and  "spikes"  towards  zero,  and  (2)  I  is  about  zero.  Visual  inspection 
of  the  figures  illustrates  this.  If  the  gain  A  is  small,  the  spike  is  hard  to  detect  but  will 
have  moor  effect  upon  the  analysis  of  the  raw  data,  since  the  first  step  in  the  analysis  is  a 
Fourier  transform,  which  is  basically  a  gain-weighted  and  phase-weighted  averaging 
process. 

The  difference  between  the  high-SNR  real  component  and  the  iow-SNR  real 
component  is  shown  in  Figure  III— 6.  The  phase  error  spikes  are  dearly  visible. 

One  of  several  possible  methods  for  locating  phase  errors  of  this  type  is  by 
inspecting  the  derivative  of  the  real  (or  imaginary)  component  with  respect  to  the  spatial 
coordinate.  If  s(i)  represents  either  the  real  or  the  imaginary  component,  then  define  the 
derivative  process 

s'(i)  =  s(i+l)-fl(i) 

which  is  a  one-sided  estimator  of  the  derivative  of  s(i).  A  spike  in  the  process  s(i)  appears 
as  a  double-sided  spike  in  s'(i)  (d.  the  engineering  "derivative"  of  the  continuous-time 
Dirac  delta  function]. 
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Figure  m-4.  Real  and  Imaginary  Pirts  (Low  SNR) 
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Figure  m-5.  Difference  Between  Real  Parts 
(High  SNR  versus  Low  SNR) 
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A  threshold  value  for  the  minimum  double-spike  amplitude  that  should  be  classified 

as  a  phase  error  can  be  derived  from  Tchebysheffs  Inequality  [14)  which  relates  the 

deviation  |X-*i|  of  a  process  X  to  the  variance  a  of  the  process  by  the  inequality 

P{|X-p|  >  ne}  <  -*«■■■ 
n^ 

irrespective  of  the  form  of  ihe  probability  density  function  of  X. 

Identify  X  with  the  derivative  process  s'(i),  so 


/i  =  E{s'} 
N 


=  ~TT  I  «'0). 

n=l 


o  =  E{(s'-p)2} 


>  Tsrl  (S'W-K)2  . 

n=l 

For  n=3,  say,  the  inequality  predicts  that  the  probability  that  |s'(i)-p|  will 
3ff  is  less  than  3~2=0.1lT. 


The  spike  detection  algorithm  is  constructed  as  follows.  Compute  the  derivative 
sequence  s'(i).  Compute  the  mean  and  variance  of  s'(i)  and  normalize  s'(i)  to  zero  mean 
and  unit  variance  by  the  transformation 
s'(i)  (- 

Choose  "n"  to  ensure  a  reasonably  low  rate  of  false  detections,  and  perform  the  test 
|s'(i)|>n  to  locate  invalid  phase  values  in  s(i).  The  test  is  performed  independently  on 
the  real  and  imaginary  components  of  s(l),  and  the  combined  results  are  used  to  select  the 
suspect  values  of  s(i)  for  which  interpolation  is  needed. 


Since  the  noisy  signal  is  oversampled  by  a  factor  of  ten  above  the  Nyquist  rate,  a 
low-pass  digital  filter  can  be  applied  to  the  measurements  to  improve  SNR. .  Averaging 


II 
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groups  of  ten  adjacent  samples  is  a  simple  approach  to  the  filtering  but  distorts  the  higher 
spatial  frequency  components  of  the  desired  signal  and  has  suboptimal  noise  suppression 
characteristics. 

A  linear-phase  filter  is  required  to  avoid  spatial  phase  distortion  of  the  measured 

signal.  A  finite  impulse  response  filter  in  the  spatial  domain  with  cutoff  at  normalized 
£ 

frequency  v  t  ~  =  -y--— —  can  be  constructed  as  follows  [151. 
cuw  'ffyquist 

We  want  a  FIE  filter  of  the  form 
m— 1 

H(z)  =  J  aj  z~*,  where 
i=0 

m  s  number  of  coefficients  (odd  integer),  and 
aj  3  filter  coefficients  (real). 

To  ensure  linear  phase,  we  define  the  filter  coefficients  symmetrically  as 
Cp  i^rO,  1, 2,  ...q, 

where 


q  =  (m-l)/2, 


c  .st. 
-i  r 


The  ideal  low-pa3s  filter  has  transfer  function  of  the  form 


V- 


1  if  v  <  v. 


cutoff 


0  if  v  >  » 


cutoff 


Taking  a  discrete  Fourier  transform  to  the  spatial  domain,  the  filter  coefficients  are 


Gain  and  phase  response  plots  (not  shown)  verify  that  this  is  a  linear  phase 


low-pass  filter  with  cutoff  near  v=0.2,  but  sidelobe  levels  arc  unacceptably  high.  Sidelobes 

are  greatly  reduced  by  using  the  smoothing  window 
0.54  +  0.46cos  (ix/q) 

0  otherwise 


wh(')  “ 


1H-16 


so  that 


CjHCJ.WjjO). 

Sidelobe  levels  of  the  filter  composed  of  the  (Figure  m — 6)  are  about  50  db  down. 
At  f=l/A,  the  highest  spatial  frequency  that  the  plane  wave  spectrum  can  contain  outside 
the  reactive  near-field  region,  the  corresponding  normalized  frequency  is  i*=0.1,  and  the 
filter's  gain  is  down  by  a  factor  of  0.92  (about  -0.7  db).  This  level  of  attenuation  will  have 
negligible  effect  upcn-the  near-field  data  analysis. 

J.  EOm^ESSOSSaLAMI^ 

Autocorrelation  (ACF)  and  power  spectral  density  (PSD)  of  the  real  and  imaginary 
components  of  the  two  scans  are  shown  in  Figures  III— 7  through  III— 10.  The  smoothing 
window  bandwidth  is  0.013^^.  No  differencing  of  the  signals  was  performed, 
although  slightly  better  PSDa  might  be  obtained.  Note  that  the  PSDs  are  down  by  20  to 
30  db  at  the  1/A  frequency.  The  noise  floor  is  down  35  db  (high-SNR  signal)  and  20  db 
(low-SNR  signal).  One  may  conclude  that  the  measured  signals  are  indeed  bandlimited  to 
1/A. 

The  composite  plots  in  Figures  D3-11  through  III— 14  illustrate  the  effect  of  the 
spike  detection,  interpolation,  and  spatial  filtering  algorithms.  The  data  traces  in  the 
figures,  from  top  to  bottom,  are 

(1)  raw  data  (real  or  imaginary), 

(2)  normalized  derivative  s' (n), 

(3)  combined  "phase  error"  flags  from  real  and  imaginary  components 
using  Tchebycheffs  Inequality  with  n=3, 

(4)  signal  after  phase  errors  are  replaced  by  interpolated  yalues, 

(5)  final,  filtered  signal. 

The  scales  of  traces  1,  4,  and  5  are  identical,  while  2  is  shown  at  0.25  sensitivity. 
Each  scan  (comprising  about  1.5  meters)  contains  1000  gain-phase  sample  pairs.  For  the 
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high  S/N  scan,  0.3  %  of  the  phase  samples  were  found  invalid;  2.2  %  of  the  low  S/N  scan 
values  were  invalid. 

PSDs  of  the  filtered  signals  are  shown  in  Figures  HI-35  and  Ht-16.  Note  that  the 
phase  error  correction  and  filtering  have  lowered  the  noise  floor  by  ten  to  fifteen  db. 
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Figure  HI-7.  Autocorrelasioa  and  rSD  of  Real  Fait  (High  SNR) 
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Figure  E-9.  Autococtiaaoc  and  PSD 
of  Real  Part  (Low  SNR) 
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CHAPTER  IV 

DESIGN  AND  CONSTRUCTION  OF  SCANNER 


A.  SCAM-SABLE 

The  scan  table  (Figure  IV-I)  has  two  movable  carriages,  each  with  one  movable 
platform.  The  two  platforms  are  referred  to  as  the  TX  and  RX  platforms.  Prior  to 
scanning,  the  TX  and  RX  platforms  start  from  known  initial  positions  defined  by  reference 
marks  on  the  scan  table.  The  initial  positions  for  the  TX  and  RX  platforms  are  near  the 
lower-left  and  upper-left  comers  of  the  available  scan  area,  respectively,  as  viewed  from 
the  operators  position.  The  stepper  motor  control  software  compensates  for  the  ofiset 
between  the  initial  TX  and  RX  positions,  so  that  coordinates  for  TX  probe,  RX  probe,  and 
test  object  are  expressed  in  terms  of  the  same  coordinate  system. 

Optoelectronic  sensors  are  located  at  the  home  positions  of  the  platforms.  When 
the  scan  table  is  initialized,  and  at  intervals  during  scanning,  the  software  moves  the 
platforms  to  the  home  positions  and  verifies  that  the  stepper  motors  are  stiil  calibrated 
(i.e.,  no  steps  have  been  lost).  If  the  software  finds  that  steps  were  lost,  the  previous 
segment  of  data  is  automatically  reacquired. 

When  near-field  antenna  measurements  are  performed,  the  antenna  under  test 
(AUT)  is  mounted  above  the  table  pointing  down,  and  the  receiving  probe  antenna  is 
mounted  on  the  RX  platform.  The  TX  platform  is  not  used  in  this  mode. 

When  near-field  bistatic  RCS  measurements  are  performed,  the  object  under  test 
(OUT)  i3  mounted,  above  the  scanning  table,  and  the  transmit  and  receive  probe  antennas 
are  mounted  on  the  TX  and  RX  platforms,  respectively. 

The  probe  antennas  are  identical  equal-length  open-ended  sections  of  X-band 
waveguide  or  pyramidal  horn  antennas.  Each  probe  is  held  in  a  mounting  bracket  that  is 
bolted  to  a  single  mounting  hole  on  its  platform.  The  TX  and  RX  probes  are  always  on  the 
platforms,  as  indicated  in  Figure  IV-1. 


Test  Object 


Figure  5V- 1 .  Scan  Table  Platforms 


Both  platforms  are  used  when  performing  bistatic  RCS  measurements;  only  the  RX 
platform  is  used  when  performing  antenna  measurements. 

C.  SQAM.PAITEBMS 

The  TX  and  RX  probe  antennas  move  through  predetermined  scan  patterns  when 
measurements  are  being  made.  The  scan  patterns  are  shown  pictorially  in  Figure  IV— 2  and 
algorithmically  in  Figures  IV— 3  and  IV-4.  During  scanning  the  probe  antennas  move 
along  the  x- axis  while  holding  y  constant. 

The  data  analysis  computer  programs  impose  some  limitations  on  the  scan  patterns. 
If  near-field  antenna  measurements  are  to  be  performed,  the  TX  scan  pattern  is 
(effectively)  a  1-by-l  scan,  and  the  RX  scan  pattern  must  be  square  with  an  even  number 
of  points  on  each  side.  If  near-field  bistatic  RCS  measurements  are  to  be  performed,  the 
TX  and  RX  scan  patterns  mu3t  be  square  and  have  an  even  number  of  points  on  each  side. 
The  x-  and  y-asmpling  intervals  for  both  the  RX  and  TX  patterns  must  be  identical. 

D.  M£ASmSMgKI-C.QmQ.URATIQM 

The  equipment  configuration  for  performing  near-field  measurements  is  shown  in 
Figure  IV-S.  The  synthesised  RF  signal  Is  taken  from  'the  output  port  of  the 
Hewlett-Packard  84G8B  Microwave  Network  Analyser  to  a  20-watt  TWT  amplifier,  and 
then  to  either  the  antenna  under  test  or  the  TX  probe  antenna  via  a  length  Lj  of 
semi-rigid  RG  .402/U  coaxial  cable.  A  length  1^  of  RG  .402/U  cable  carries  the  received 
signal  from  the  RX  probe  to  a  preamplifier  and  then  to  the  input  port  of  the  84038. 
There  is  additional  amplification/attenuatiou  in  the  test  channel  gain  G,c  setting  in  the 
8408B,  which  may  be  adjusted  differently  on  different  measurement  runs.  Lj  and  are 
typically  both  25  feet  long. 

Gaia  and  phase  are  calculated  continuously  by  the  S408B  and  are  captured  by  the 
dual— cfcannei  A/D  converter  for  the  IBM  PC.  The  analog  processing  circuits  in  the  840SB 
can  operate  at  either  100  hz  or  10  khz  bandwidth.  The  10  khz  bandwidth  is  used  to  avoid 


smearing  the  measured  data  during  fast  scanning. 
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Scan  Table 

coordinate  system  origin 
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c 
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c 
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20 


Pseudo-code  for  Antenna  Measurement  Scanning 


NXRJCNYRX  “  number  of  points  along  x,y  for  RX  probe 
ECRXJYRX  ■>  loop  indices  for  RX  scanning 
XRXORG.YRXORG  *  lower  left  comer  of  RX  scan  pattern 
XRX.YRX  m  current  coordinates  of  RX  probe 
IRXO  »  RX  probe  orientation 
DELTXY  »  sample  interval  (meters) 

DO  20 IRXOU 
DO  lOIYRX-lJfYRX 
DO10IXRX-1JNXRX 
XRX»XRXORG+(IXRX— 1  )*DELTXY 
YRX«=YRXORG+(IYRX-l)‘DELTXY 
CALL  MOVE(XRX,YRX) 

CALL  MEASURE(GAIN,PHASE) 

CALL  WRTTH(GAIN>PHASE) 

CONTINUE 

PAUSE  'Change  probe  orientation  now.’ 

CONTINUE 
STOP  Done.' 

END 


F.gure  FV-3.  FORTRAN  Cede  for  Antenna 
Measurement  Scan  Pattern 
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C  Pseudo-code  for  Bistatic  Measurement  Scanning 
C 

C  NXTXJJYTX  *  number  of  points  along  x,y  for  TX  probe 
C  IXTXJYTX  m  loop  indices  for  TX  scanning 
C  XTXORG.YTXORG  *  lower  left  comer  of  TX  scan  pattern 

C  XYX.YIX-  current  coordinates  of  TX  probe 

C  rrxo  *  TX  probe  odentarion 

C  NXRX^TYRX  *  number  of  points  along  x,y  for  RX  probe 
C  DCRXJYRX  -  loop  indices  for  RX  scanning 
C  XRXORG.YRXORG  *  lower  left  comer  of  RX  scan  pattern 
C  XRX.YRX  *  current  coordinates  of  RX  probe 
C  IRXO  «  RX  probe  orientation 
C  DELTXY  »  sample  interval  (meters) 

C 

do  30  rrxo«u 

DO  30 IRXO-U 

DO20IYTX-1.NYTX 

do  20  ixtx«i,nxtx 

XIX-XIXORG+(DCrX-l)*DELTXY 

rrX-YTXORGv{IYTX-l)»DELTXy 

CALL  TXMOVECXTXYTX) 

DO10IYRX-1JSTYRX 
DO  10  DfRX-1  JOCRX 
XRX-XRXORG+<IXRX-l)*DELTXY 
YRX««YRXORG+(IYTX-l)*DELTXy 
CALL  MOVECXRX.YRX) 

CALL  MEASUKE(GAIN,PHASE) 

CALL  WRITE(G  AIN  .PHASE) 

10  CONTINUE  - 
20  CONTINUE 

PAUSE  'Change  probe  orientations  now.' 

30  CONTINUE 
STOP  Done.' 

END 

Figure  IV-4.  FORTRAN  Code  for 
Bistatic  Scan  Pattern 
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E.  DATA  ACQUISITION  SYSTEM 

Data  acquisition  is  controlled  by  a  computer  program  that  runs  on  the  IBM  PC. 
The  program  controls  the  scanning  table  motors  and  A/D  converter  and  stores  the 
measured  data  to  disk. 


The  conversion  time  of  the  12-bit  A/D  converter  is  typically  25  /s.  Gain  and  phase 

are  sampled  simultaneously  by  dual  sample/hold  amplifiers  and  then  digitised  in  succession 

by  the  A/D  converter.  The  measurements  are  the  gain  (V  )  and  phase  (V £  voltages  from 

the  network  analyzer.  The  actual  gain  and  phase  are  computed  as  G  =  V  -f  .  and 

g  gam 

where  fgain  is  a  gain  calibration  factor  (20.0  db/volt)  and  is  a  phase 
calibration  factor  (100.0  degrees/volt)  for  the  840SB  analyzer.  The  gain  and  phase  are 
combined  into  a  complex  number  and  stored  in  a  disk  file. 


CHAPTER  V 

COMPARISON  OF  MEASUREMENTS  WITH  SOLUTIONS 


A.  INTRODUCTION 

The  far-field  RCS  that  is  calculated  via  near-field  RCS  measurements  may  be 
checked  by  comparing  it  to  the  far-field  RCS  predicted  by  electromagnetic  theory.  The 
conducting  sphere  and  disk  are  convenient  objects  for  this  comparison,  since  an  analytic 
solution  exists  for  the  sphere  and  a  physical  optics  solution  may  be  applied  to  the  disk. 
Equations  for  scattering  from  sphere  and  disk  may  be  found  in  Appendices  C  and  D, 
respectively.  In  this  chapter  we  present  a  quantitative  comparison  of  the  far-fieid  RCS 
obtained  via  near-field  measurements  with  the  calculated  far-field  RCS  of  the  conducting 
sphere. 

Although  full  bistatic  scattering  measurements  were  made  on  the  sphere  and  the 
disk,  funding  and  time  limitations  precluded  our  performing  the  data  reduction  and 
analysis  of  that  data.  The  results  presented  here  represent  scattering  from  the  sphere  only. 

These  near-field  measurements  were  made  at  10  Ghz  using  the  experimental  setup 
described  in  Chapter  III.  The  target  is  a  precision  6  inch  diameter  aluminum  sphere 
mounted  above  the  scanning  table.  The  measurements  were  made  on  a  64  by  64  grid  of 
points  with  the  transmit  antenna  stationary  and  directly  under  the  target,  so  that  the 
target  illumination  was  essentially 

E=|E|ey. 

B.  BACKSSQMP..£.UBlRAgJIPJ 

In  figure  V— 1  may  be  seen  the  magnitude  of  the  co-polarization  component  of  the 
raw  data  measured  by  the  receiving  antenna.  The  target  sphere,  located  at  (x,y) 
coordinate  (1.02,1.02)  meters,  is  above  the  lower  right  edge  of  the  plot,  corresponding  to 
the  peak  amplitude  of  the  raw  signal. 

Figure  V-2  shows  a  background  plot  with  all  parameters  identical  to  the  previous 
plot  with  the  exception  of  the  target  sphere,  which  was  removed  from  the  test  volume. 
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Figure  V-2.  Co-polaricadoa  Background  Magnitude 
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Although  there  is  no  coherent  background  signature  visible  In  the  figure,  the  raw  data  are 
noticeably  cleaner  after  this  background  scan  is  subtracted  coherently  from  it  (Figure 
V— 3).  The  improvement  in  signal/noise  ratio  is  even  more  apparent  in  Figures  V— 4  and 
V-5,  which  are  contour  plots  of  the  imaginary  component  of  the  co-polarization  before  and 
after  coherent  background  subtraction,  respectively.  Similar  improvement  occurs  in  the 
cross-polarization  components. 

Coherent  background  subtraction  is  performed  for  all  data  shown  in  the  figures  in 
the  remainder  of  this  section. 


The  first  step  in  the  analysis  of  the  raw  data  is  a  two-dimensional  discrete  Fourier 
transform  (DFT).  The  signal  processing  aspects  of  the  data  analysis  are  illustrated  by  the 
effects  of  several  tapering  methods  that  were  tried.  Figure  V— 6  shows  the  magnitude  of  the 
DFT  of  the  co-polarization  signal  with  no  tapering;  this  is  equivalent  to  a  rectangular 
("boxcar")  taper  in  the  continuous  domain.  The  direct  scattering  from  the  target  shows  up 
as  the  large  peak  in  the  figure,  and  probe-to-probe  coupling  appears  at  the  upper  left  and 
lower  right  edges  of  the  figure.  The  side-lobe  level  due  to  leakage  from  the  main  peak  is 
about  30  db  below  the  peak. 

In  figure  V-7  is  shown  the  magnitude  of  the  DFT  of  the  co-polarization  signal 
after  a  separable  cosine  taper  of  the  form  1  +2C03^  ~  was  applied  to  both  the  x-axis 
and  the  y-axis.  Note  that  the  sidelobe  level  is  better  than  40  db  down  and  the 
probe-to-probe  coupling  peak  is  much  more  localized,  although  there  is  little  fine  structure 
discernible  on  the  peaks. 

figure  ¥-8  shows  the  magnitude  of  the  DFT  of  the  co-polarization  signal  after  a 
separable,  partial-cosine  taper  was  applied  to  both  axes.  This  taper  consists  of  a  cosine 
taper  spliced  onto  a  boxcar-type  taper,  so  that  the  middle  half  of  the  data  are  unaffected 
while  the  first  and  last  quarters  of  the  data  set  are  tapered.  The  sidelobe  suppression  ol 
this  taper  is  generally  poorer  than  that  of  the  cosine,  but  for  this  data  set  is 
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Figure  V-6.  Co-poltrratioa  DFT  Without  Taper 
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comparable  to  that  produced  by  the  cosine  taper.  Much  more  structure  is  visible  on  the 
peaks  (a  common  characteristic  of  this  type  of  taper). 

A  similar  judgment  may  be  made  based  on  the  calculated  scattering  matrix. 
Figures  V— 9  through  V— 11  show  the  magnitude  of  calculated  using  no  taper,  a  cosine 
taper,  and  a  partial  cosine  taper,  respectively.  The  cosine  taper  has  the  lowest  sidelobe 
levels,  while  the  partial  cosine  taper  has  higher  sidelobes  but  more  derail  in  the 
transformed  data.  In  general,  the  additional  detail  makes  the  partial  cosine  taper 
preferable  to  the  cosine  taper. 

D-  PROBE-PROBE  COUPLING 

Figures  V— 3  and  V— 8  illustrate  the  manner  in  which  the  measured  data  are 
transformed  from  a  Cartesian  coordinate  system  with  units  of  meters  (Figure  V-3)  to  an 
angular  spectrum  with  units  of  reciprocal  meters  (Figure  V-8)  in  k-space.  The  directly 
coupied  signal  from  the  transmit  probe  to  the  receive  probe  appears  in  the  angular 
spectrum  as  a  broad  peak  with  incident  wavenumbers  near  the  horizon  (i.e.,  |k  j  g  0)  in 
angular  space. 

This  direct  signal  can  be  separated  from  the  desired  signal  if  the  desired  signal  is  not 
too  close  to  the  direct  signal  in  k-space.  Note  that  (Figure  V-8)  the  probe-probe  coupling 
appears  as  a  broad  peak  at  the  top  left  edge  of  the  plot  and  spills  over  (due  to  wraparound 
in  the  DF1)  to  the  lower  right  edge.  This  broad  peak  is  distinct  from,  and  do®  not 
corrupt,  the  desired  signal  in  these  plots. 

Under  these  conditions,  the  only  detrimsitsl  effect  of  probe-probe  coupling  is  that 
it  increases  the  required  dynamic  range  of  the  receiver. 

E.  REGION  OF  DEFINITION 

The  region  of  definition  of  the  reconstructed  RCS  is  limited  by  the  size  of  the  scan 
area  relative  to  the  apparent  angular  extent  of  the  object  as  viewed  from  the  center  of  the 
measurement  scan.  In  general,  the  actual  area  of  definition  is  relatively  small  when  viewed 
in  k-space.  The  reconstructed  phase  of  Sg  shown  in  Figure  V-12  is  valid  only  over  a  small 
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semicircular  region  that  is  approximately 
-50  meter-1  <  kx  <  50  meter-1 
60  meter-1  <  ky  <  140  meter-1. 

The  actual  region  of  validity  is  blurred  due  to  the  sampling  geometry  and  signal  "leakage" 
between  adjacent  DFT  values  The  leakage  has  been  reduced  in  these  plots  by  applying  a 
smoothing  taper  before  performing  the  Fourier  transforms.  The  region  ox  validity  appears 
in  the  figure  as  an  area  of  smoothly-varying  phase  that  is  distinct  from  the  phase  noise 
that  covers  the  remainder  of  the  k^-ky  plane. 

F.  CSIXIMS 

The  region  of  definition  appears  also  in  cut  plots  through  the  k  ~k  plane.  A  cut 
through  Figure  Y— 12  along  the  line  kx=0  yields  the  plot  shown  in  Figure  Y-  *13.  The  phase 
function  is  well-behaved  for  60  meter-1  <  ky  <  140  meter-1,  while  it  oscillates  rapidly 
and  randomly  outside  that  region. 

Figures  V— 14  through  V-26  show  overlay  plots  of  the  predicted  and  measured 
Sx,  and  Sy  in  cuts  along  the  kx  and  ky  axes.  A  fixed  offset  of  6  db  has  been  added  to 
the  measured  scattering  magnitude  to  facilitate  comparison  of  predictions  and  measured 
data  in  these  plots.  This  offset,  occurring  consistently  re  all  of  the  measured  scattering 
data,  is  attributed  to  a  calibration  inconsistency  that  time  limitations  prevent  our 
resolving  Recall  that  the  region  of  definition  in  these  plots  is  approximately 
-50  meter-1  <  kx  <  50  meter-1 
60  meter-1  <  ky  <  140  meter-1. 

Figure  V— 14  show?  the  predicted  and  measured  magnitude  of  the  function  Sff  along 
a  cut  kx  *  0  meter-1.  The  deep  notch  in  the  predicted  magnitude  is  a  point  of  phase 
ambiguity  in  this  coordinate  system  representation.  Over  the  region  of  definition,  the 
reconstructed  magnitude  is  essentially  flat  with  about  ±2  db  of  residual  ripple. 

Figure  V— 15  shows  the  predicted  and  measured  magnitude  of  the  function  along 
—1 

a  cut  k  =  0  meter  *  Along  this  cut,  the  predicted  S..  is  zero  except  {or  an  impulse  at  fhe 
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origin.  The  reconstructed  funeticn  is  about  -20  db  from  indicating  that 
cross-polarisation  leaksge  occurred  due  to  antenna  misalignment. 

Figure  V-J.6  shows  the  predicted  and  measured  magnitude  of  the  function  Sy  along 
a  cut  =  0  moter-1.  'Che  predicted  Sy  is  accurately  reconstructed  with  about  *2  db  of 
ripple. 

Figure  V— 17  shows  the  predicted  and  measured  phase  of  the  function  S g  along  a  cut 
kx  =  C  meter-1.  In  these  phase  cut  plots,  it  is  important  to  view  only  the  actual  region  of 
definition  gWen  above.  Over  that  region,  the  phase  agreement  >s  reasonably  good;  the 
upward  parabolic  slope  to  the  right  is  discernible. 

Figure  V-18  shows  the  predicted  and  measured  magnitude  of  the  function  5y  along 
a  cut  kx  =  0  meter- ^ .  Again,  over  the  region  of  definition,  the  reconstructed  phase  agrees 
reasonably  veil  with  *he  predicted  phase,  provided  that  one  takes  into  account  the  phase 
wrap  that  occurs  at  about  U5  meter-1. 

Figure  V-t9  shows  the  predicted  and  measured  magnitude  of  the  function  Sp  e'ong 
a  cut  ky  =  65.1  meter-1.  Along  this  cut,  where  the  region  of  definition  is  approximately 
Iky!  i  50  meter-1,  the  agreement  between  the  two  plots  is  excellent,  with  about  2  db 
maximum  error. 

Figure  V— 20  shows  the  predicted  arid  measured  magnitude  of  the  function  along 
a  cut  kx  =  65.4  meter-1.  Agreement  is  again  quite  good,  with  the  same  offset  of  6  db. 
Tbe  same  may  te  said  of  Figures  V— 21  and  V-22  whwh  show  the  predicted  and  measured 
magnitude  of  the  functions  5  and  S  along  the  same  cut. 

Figures  V— 25  through  V-26  show  tbe  predicted  and  measured  phases  of  the 
functions  S^,  Sx,  and  Sy  along  a  cut  =  65.4  roster-1.  The  measured  phase  plots 
show  a  linear  additive  phass  term  with  slope  of  about  0.3  degree/meter,  this  term  is 
probably  due  to  target  misalignment.  There  is  in  addition  a  phase  inversion  of  S  .  and 

Y 

Sx which  affects  the  maimer  in  which  those  two  functions  change  as  the  origin  (kx  =  0)  U 


C-.  SUMMARY 

In  this  section  we  presented  both  predictions  and  near-field  measurements  of 
tar-field  bistatic  scattering  of  the  conducting  sphere.  The  hey  points  of  the  presentation 
may  be  summarized  as  follows: 

(L)  The  magnitude  of  the  measured  scattering  agrees  well  with  the  predicted 
scattering  within  the  region  of  definition  of  the  k-space  reconstruction.  The 
discrepancy  between  measurement  and  prediction  consists  of  a  ±2  db  peak 
ripple  superimposed  upon  a  constant  calibration  offset  of  6  db. 

(2)  Agreement  of  the  scattering  phase  is  not  as  good.  This  is  due  primarily  to 
two  effects: 

fa)  Phase  measurements  are  particularly  sensitive  to  antenna  axial 
alignment  perpendicular  to  the  plane  of  scanning. 

(b)  "Wrap-around"  occurs  in  the  phase  of  the  measured  data  due  to 
(fixed)  phase  offsets  in  the  network  analyzer. 

(3)  Coherent  background  subtraction  is  an  effective  technique  for  improving 
signal— to-coise  ratio. 

(4)  Spatial  filtering  to  remove  unwanted  scattering  from  the  envmonmeat  may 
be  performed  in  the  k-space  Fourier  transform  representation  of  the 
measured  data. 

(5)  Sidelobes  iu  the  Fourier-transformed  measurements  are  greatly  reduced  by 
smoothing  the  measurements  with  a  smoothing  window  prior  to  per  Vicing 
the  Fourier  transforms. 

(6)  Probe-to-probe  coupling  is  readily  separated  from  the  desired  target 
scattering  >cform<ttion  in  the  k-space  representation  of  the  measured  data. 


CHAPTER  VI 

DATA  REDUCTION  INVESTIGATION 


A.  INTRODUCTION 

The  existing  algorithms  for  computing  near-field  antenna  patterns  and  radar  cross 
sections  of  objects  impose  significant  computational  burdens.  In  this  section  we  describe 
modifications  and/or  alternatives  to  the  existing  algorithms  which  offer  the  potential  for 
reducing  the  number  of  computational  steps  needed  to  obtain  useful  near-field 
measurement  data.  The  investigations  concentrate  on  three  primary  computational 
alternatives  described  in  the  Mowing  sections.  The  following  other  computational 
alternatives  are  not  considered  here;  (1)  adaptive  nonuniform  sampling,  (2) 
non-rectangular  lattice  (or  grid)  sample  spacing,  (3)  statistically  reduced  number  of 
near-field  data  used  in  ti>e  far-field  transform,  and  (4)  derivative  data  measured  in  the 
near-field. 

Two  approaches  for  reducing  the  burden  of  the  near-field  computations  are 
considered  here.  The  first  approach  io  based  on  writing  the  equations  for  computing  the 
scattering  matrix  directly  in  "ector  and  matrix  notation.  The  current  equations  cumpute 
e  scattering  matrix  in  several  steps.  Each  step  requires  computations  and  transfers  to 
data  storage  which  impose  computational  burdens.  These  burdens  can  be  reduced  by  using 
a  direct  matrix  formulation  cf  the  scattering  equation  In  addition,  the  other  alternatives 
for  reducing  either  the  amount  of  measured  near-field  data  or  the  computations  necessary 
to  determine  the  scattering  matrix  are  more  convenient  to  work  with  in  terms  of  the  direct 
matrix  formulations. 

The  tecond  approach  '•«  based  on  expansions  of  the  near-field  functional  components 
which  allow  us  to  impose  constraints  on  the  near-field  functions  The  constraints  can  be 
based  on  (1)  the  desired  or  expected  analytical  behavior  of  the  functions,  (2)  the  physics  of 
i‘!v  measurement  environment,  or  (3)  the  information  content  of  the  data  process. 


B. 


In  investigating  the  computations  from  a  matrix  perspective,  we  first  rewrite  the 
equations  for  the  scattering  matrix  in  vector  matrix  notation.  Expressed  in  this  form,  the 
scattering  equation  is  more  easily  interpreted  and  adapted  to  improve  the  efficiency  of 
computing  the  relationships  between  the  scattering  components.  Next  we  relate  the 
scattering  matrix  formulation  with  the  vector  formulation  of  the  first  subsection  and 
compare  the  results. 

The  explicit  solution  for  the  scattering  matrix  Xjj(K,L)  in  terms  of  the  measured 
data  array  parallels  the  solution  given  in  Chapter  II  except  that  this  derivation  solves  for 
the  scattering  matrix  in  terms  of  the  measured  data  vector  and  matrix  representations. 
The  equation  development  is  descrioed  in  terms  of  the  continuous  Founsr  transform  for 
exposition  convenience.  Equivalent  results  can  be  obtained  fo>  the  point  measurement 
computation  of  the  Discrete  Fourier  transform  using  the  Fast  curier  transform. 

The  bistatic  measurement  equation  is 

•f«  ,  , 

"dkdl 

'  ‘  '  JJ . 

where 


i(I)  •  2 


J‘x(k)  =  the  vector  transpose  of 
k  =  scattering  (emerging)  direction, 

1  =  incident  (illumination)  direction, 


f  j  =  +  Pyj®y  +  —  receiver  probe  coordinates, 

?2  =  Pxtex  +  Pyt®y  +  Pztez  =  transmit  probe  coordinates. 

In  the  laboratory  we  measure  a  multidimensional  scalar  complex  array 
M(Pxr>Pyr-Pxt-P  yj.fj.fj)  at  discrete  points  whose  components  elements  are  the  measured 
values  of  b0(fj,f2) 

(Pxr>Pyr'Pzr)  6  {^xr^yr^zPj-  Pzr  =  COflStant- 
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b0(fi,?2)e  1  2  df !  df2  * 


(PxfPyt-Pzt)  e  {(Vyt-^t)}*  Pzt  =  constant’ 

=  {it  2}  =  an  index  denoting  receive  probe  orientation,  and 
{  =  {1, 2}  s  an  index  denoting  transmit  probe  orientation. 

In  the  first  step  of  the  solution,  we  take  the  inverse  Fourier  transform  with  respect 
to  the  respective  transmitter  and  receiver  locations  f  j,  fj  of  both  sides  of  the  equation. 

i‘i(k)-T11(k,i)-i10(i)' 

e  2  1  didl  e  1  zdf1df2 

=  (2x)6  (£)  -  TX1  (£, I)  •  1 10(I)  -  «Cm+l)  -  <5(5-k)  dk  dl 

=  (2 *)6 

where  k  is  defined  in  terms  of  its  wavenumber  k^,  kyj,  ka  components  in  the  coordinate 
system  of  Tu(ic,l).  Note  that  (t  and  ^  are  constants  in  the  above  equation.  We  now 
write  the  equation  in  terms  of  the  k,I  arguments 
B(k,-I)  =  (2r)6  ^(k)  -Tn(J,I)  -!10(I) 

where 


-  v  ff  -ik-fj+iT-f,  ,  , 

B(k,-I)=  b0(f1,f2)*e  2 

The  minus  sign  of  the  vector  1  indicates  the  direction  of  the  incident  radiation 
relative  to  the  scattering  object.  In  our  laboratory  measurements  the  elements  of  M(-)  axe 
built  up  by  a  concatenation  of  scans  by  the  rec e  ve  probe;  during  each  scan  we  vary 
and  P  and  hold  ?u,  Pxt,  Pyj,  £r,  and  constant.  If 
«  number  of  discrete  value  of  P^, 

N^  =s  number  of  discrete  value  of  P^, 

=  number  of  discrete  value  of  P^, 

Ny.  =  number  of  discrete  value  of  Pyt, 


then  M(-)  consists  of  4NxtN^  matrices,  each  of  which  contains  N^N,^  complex  scalars. 
The  factor  of  4  comes  from  (r=  {1,2}  and  £t=  {1,2}. 

The  vectors  i01(£)  411(1  ^io$)  ^ve  s?e<^c  components  reflecting  their 
dependence  on  the  probe  orientations  ((r=  {1,2},  {1,2})  polarisation  components  (5,0) 

of  the  transmitted  and  received  signals.  Next  we  define  two  2x2  matrices  IK  and  IE 
composed  of  the  respective  ig^(ic)  and  1^(1)  vectors  in  terms  of  their  specific  function 
arguments  £,  0  and  0,  thus 

y£_  'oi^’^rl)  _  ^l^’^rl^  *O10&4l) 

^01  ■  !01^’^  !O10^’^  • 

^l^’W  _  ^O^’^tl^  IlO0^*’^tl^ 

■  W  I1O0^’^ 

The  quantity  B(k, — I)  is  also  a  function  of  the  probe  and  polarization  arguments  so 
we  define  5,  a  2x2  matrix  with  the  diagonal  terms  specified  by  the  respective  values  of 
B’(-)  for  each  of  the  four  vector  combinations  of  the  IK  (IE)  matrices.  Using  the  above 
definitions  we  can  write  the  scattering  matrix  equation  as 

K(k,I)  -  (2x)6IK(k)-T11(k,I).n:(I)t. 

If  the  matrices  IK  and  IE  have  rank  two,  their  inverses  exist  and  we  can  write  the 
scattering  matrix  as 

Tn(iU)=(IK)-I.K.(IEtr1/(2x)6. 

The  matrix  inverses  for  the  IK  and  E?  matrices  are 

(TK)-l  =  "W^rl)  1/A(r) 

r01^k'4l) 

IO10^1’^ 

_IQ10^^tl^  ^l/'^tP 
where  A(r)  and  A(t)  are  respectively 

^(r)  =  ^l^’^rP  ^010^’^  ~*O10'k,^rP  *0l/k,'riP 

=  ho^’^tl)  I1O0^1’^  ”  IlO0^,'tl^ 


The  E  matrix  is  given  by 

2  _  Brl  tl  Brl  t2 
■  Br2  tl  Br2  t2 

The  subscripts  ri  and  tj  (i,j=l,2)  in  the  Brf  tJ.  terms  identify  the  orientation  of  the  fy 
probe  components.  We  can  then  write  the  scattering  matrix,  T^j(k,i),  explicitly  in  terms 
of  the  vector  elements  and  scaling  components,  A(r)A(t)(2x)  ,  as 


tl~I01<^’{rl)Br2  tl  W&^rt^rl  t2  W&WBr2  t2 
-IQlo(k,fr2)Brl  tl+I01  fl^rl^r?  tl  ^01  d&^ri^rl  t2+I01  ^’^rPBr2  t2 


*01  “*01^1'W 


(A(r)A(t)(2x)6)-1 


We  can  now  write  the  elements,  Tjj(id=l,2)  of  the  Tjj(k,I)  matrix  as 

TU  =  Poi^’^rt^rl  tl-I01^^ri^r2  tl^Ol^’W  ~ 

t2-W^’*rl)Br2  t2%^7^tlWA^A-t^2^  ^ 


T12  _  ^-I01^^’^Brl  tl-I0U^rl)Br2  tl^Ol^’^P  + 

^01^^’^Brl  t2-I01^^rl^Br2t2^I01^’^tl^^AWA^t^2^  ^ 


T21  =  (-I01^^^r2^Brl  tl+I01^’^rl^Br2  tl^Ol^’^  “ 

(-I01fl^’WBrl  t2+!O10^'ri)Br2  t2%^*’^tlWA(r)A^^2^  ! 
T22  =~^-I01(/^’'r2^Brl  tl+I01^’’rl^Br2  tl^Ol/*’^  + 

(-W^r2>Brl  t2+I01^rl)Br2  t2)I0l/7^tl)KA(r)A(t)(2x)6)-1. 


The  T-j  components  of  Tj^k,!),  for  particular  values  of  (^,^,1,^,1^),  are  the 
transfer  function  that  maps  incident  $,<p  polarizations  into  scattered  (emerging)  0,<j> 


polarizations. 

We  now  show  the  equivalence  of  the  scattering  matrix  as  derived  above  and  the 
matrix  representation  derived  in  Chapter  U.  We  Srst  derive  the  explicit  solution  for  the 
scattering  matrix  T^fKjL)  in  terms  of  the  measured  data  array.  The  bistatic 
measurement  equation  is 
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D&fj) 


where  the  equation  elements  are  as  defined  in  the  preceding  sections  and  D(k,r^)  is  a  scalar 
coupling  product. 

In  the  laboratory  we  measure  a  multidimensional  scalar  complex  array 

VPxr-Pyr’Pzr,Pxt-PyfPZt^r> «t>  at  di3Crete  P°ints 

(Pxr- W  e  {^xr^yr^ar)}-  Pzr  =  C0“toBt’ 

(Pxt>Pyt-Pst)  6  {(^ytjM*  Pzt  =  «“> 

=  {1,  2}  =  an  index  denoting  receive  probe  orientation,  and 
=  {1,  2}  =  an  index  denoting  transmit  probe  orientation. 

In  the  first  step  of  the  solution,  the  scalar  coupling  product  D(-)  is  computed  as  a 
phase-shifted  two-dimensional  Fourier  transform  of  B(  ■ )  with  respect  to  the  transmitter 


coordinate  r2: 


— im*r<.  .  ,  — il*r.  4  , 

b0(r1,?2)  e  "dr2=  4(k).  Tn(k,l)-T10(l)  e  dk  dl 


-il-r.  .  „  i(k-m)-r2  „ 

i  *  ill  a 


which  we  write  as 


,4  4 

©(m,^)  =  4(5)  •  J  Tu(k,i)  -i10(I)  e  11  fl  dk  dl  (2*)2  *(k-m) 
n(5,?i)  =  (2x)2  4(5)  j  tu(5,I)  -i10(I)  e*11  ri  dl. 

C(S,r  j)  is  now  written  as  a  matrix  to  reflect  its  multi-component  structure. 

In  our  laboratory  measurements  the  data  components  l>0(rpr2)  toiHf  up  as  a 

concatenation  of  scans  by  the  receive  probe;  during  each  scan  we  vary  and  Pyr  and 

hold  P  ,  P,„  P  rt,  £  ,  and  constant.  Each  U(k,r.)  component  consists  of  four  complex 
zt  xt  yi  *r  i  * 

scalar  corresponding  to  the  combinations  of  fr={l,2}  and  <ft={l,2}. 
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In  the  second  step  of  the  solution,  the  scattered  spectrum 


is  computed  from  U(k,  r^)  and  the  receive  probe  response  function  Tqj(-)  by  means  of  a 
simultaneous  linear  equation  solution.  From  a  prior  equation  is 

U(^;i)  =  (2x)2T‘1(k).F(k,?1). 

Expanding  explicitly, 

D(k,rj,fr,fj)  =  lgi^(k,fr)  +  I(ji^(k,{r)  ^^(k,rj,fj). 

To  solve  for  F(-)  we  need  two  equations  obtained  by  measuring  D(-)  at  two  different 
D(k.ry^rli?t)  =  W^)  Fj(k,rp4{)  +  lQi^(k,{rl)  F^(k,rj,ft) 

^(k.rj,^,^)  =  ^01  "*■ 

Then 

F5(k,rj,{t)  =  ^(k^.^.f^Ig^k,^)  ~D(k,r^,fr2,ft)lQj^(k,{rj)jA(r) 

^(k.rp^)  =  [n(k>r1,4rl,ft)Ig1^(k1^2)  -D(k,r1,ft2>ft)lQ1^(k,frl)JA(r) 
with  non-zero  determinant 
A(r)  s 

For  each  transmit  probe  position  (P^.P^.P^)  and  orientation  ft,  F(-)  is  defined  over  all 
valid  k.  Note  that  F(-)  is  not  a  function  of  {f.  A  particular  realization  F(-)  is  the 
scattered  spectrum  from  the  object,  given  a  particular  position  and  orientation  ft  of  the 
transmit  probe. 

In  the  third  step  the  vector  coupling  product  is  calculated  as  a  two-dimensional 
phase-shifted  Fourier  transform  with  respect  to  the  transmit  probe  coordinates.  From  a 
prior  equation 


we  define 
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S&D-Tn&lj.ljotf). 


and  write 


-il*r 

P(5,?l)=  TTtf,I)-e‘a  'ridl. 


The  equation  is  now  Fourier  transformed  with  respect  to  r^ 

*  „  -in-r-  „  rr  -ifl+nj-r,  „  „ 

F(k,r1)-e  1  drj  =  £J(k,l)-e  1  dl  drp 

F(U)  =  (2x)2  S(U)-dI  i(!+n) 

F(k,I)  =  (2»)2Q(k,-i) 

This  vector  coupling  product  represents  the  8,$  components  of  the  scattered  spectrum  as  a 
function  of  incident  I,  emerging  k,  and  transmit  probe  orientation  ft. 

The  fourth  and  final  step  is  to  calculate  the  scattering  matrix  Tn(M).  The 
equation 

q&r1)-Tu(ilI)-l10(I) 

can  be  written  explicitly  in  terms  of  the  respective  ^terms  as 

Q^(k,ri,fti)  Qg(k,rp^t2)  _ 

■  Q^(k,rp5{2) . 

[  Tn^l.^tl)  Tll^^’ft2)  1 


L  Tll^ri»?tl)  Tll,A2>  -I  i  ho^^’^tl) 

Solving  for  Tjj(k,I)  we  write 

T11(i,i)=^(k,?pet).T10(i,{tr1 

where 

t  t  v1  -  -*108(^2)  AM-1 

w  1%ftw  ^.wr' 

with  non-zero  determinant 

M*) 5  hoffiAilhotfaa)  ~ 

which  results  in 
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QflO^pCti)  Q^>ri>^2^ 


A(t) 


-1 


-I 

We  can  now  write  the  components  of  5^(2,!)  as 

^  /it  T\  _  »ft2>  ’ 

Tllflfl(*’1' - — -  A(t) - C - 

Tii^&i) =  ^i^hs£hk 

t,  „*t»  _~  V^^ti^io^’^  + 

An  <j>e^lj  - 

^  ^  _  ~  q^(^,eti)ii0tf<T  ,^t2) + 


t-nw~''* - sit] 

The  four  components  of  Tjj(k,l),  for  particular  values  of  (2,1),  are  the  transfer  function 
that  maps  incident  9,<f>  polarizations  into  scattered  (emerging)  9,<f>  polarizations. 

To  compare  the  vector  and  above  formulations  of  the  scattering  equations,  it  is 
necessary  to  substitute  for  the  respective  $,  F  and  U  matrices  io  show  the  equation  in 
terms  cf  the  measured  B'  data  set. 

We  have,  from  a  preceding  section, 

F(2,l)  =  (2x)2$(2,-1). 


We  use  a  previously  defined  equation  relating  U  and  F  to  write  the  T  components  in  terms 
of  the  IJ  and  I  matrices,  i.e., 

n(S,;1)=(2x)2T*1(5).F(i;,;1). 

The  Fourier  transform  with  respect  to  r^  and  the  wavenumber  argument  1  is 

B(k,l)=(2r)2T‘1(2).F(i,l). 

Thus,  we  obtain 

D(k,I)=((2T)4T‘1(k)).U(k,4). 

Recall  that  EJr^  r2)  is  defined  in  terms  of  C(k,  ?j)  as 


Taking  the  Fourier  transform  of  both  sides  with  respect  to  r^  and  ?2  we  have  . 


„  -lm-rj-in-r,  ,  « 

yrPr2) e  dridr2’ 


_  ,  ,  -ik-r9-im-r,-in*r9  „  „  „ 
UCkjTj)  e  2  1  2dkdr1dr2. 


The  results  of  the  transform  operations  yield 

* 

“ira*r. 


H(m,n)  - 


17(5,^)  *e  1  dk  drj*(2r)  ‘$(k+n) 


or 


S0(k,l)*(2  t)2.U(M). 

Using  the  equation  for  TJ  determined  earlier,  we  now  have 

s(U)=(2x)6.4(£)).3(U). 

We  can  now  write  the  components  of  in  terms  of  the  E  components  as 

U(£,T)  -lJ1(Jr1.H(S,l)/(^e. 

The  I  and  E  elements  of  IJ  are 


4ft”1- 


A(r)~ 


!oi^k’^  -WHi) 

where  A(r)  is 


and 


Bllu  Efit2 


u  u 
*•  Br2tl 


r2 12 


We  can  now  write  {?(£,¥)  as 

U(k  I)  =  f  <^k,ri,^ti) 
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(Wk’^Rn  tl  "WMA  tl'  (l01^k’^Brl  t2  -I01^k,^rl)Br2 t2^ 
Hqi  ‘r2^Brl  tl  J'!31  0&£rl-Br2  tP  Hoi 5^^Brl  t2  +I01  ^k’^rPBr2  t2^ 


x((2r)6A(r))-1. 


We  can  now  write  the  components  of  (it, ! )  aa 

™  ,?1s 

T110^k’1)  "  . A'(VJ 

Tll^W) = — - z — - - 


— »  / 

ln<pe(Kl> =  Ipr 


LiW 


,(k,i) 


—  v-/ 

-  ^*.kT,ftl)I10j(I,(t2)  +  ^(kJifg^io^l’ftl) 

i»)  5  "•  ■  A  If  1  • 


W 


When  the  vclues  for  tl  s  fj  elements  are  substi  uied  into  the  equations  for  the  ¥  elements 
we  see  that  th.  scattering  matrix  yalues  are  identical  to  those  derived  earlier.  Thus,  the 
representation  -re  will  use  for  computing  the  scattering  matrix  is 

TU(U)  -  (IKtf; r1  •  B(k , J) . (tc (1 ) lrV(2r)6 

where 

(TKT1  =  1/A(r) 

'-I0l/k,W  ;pid(k-fri)  J 


(TE*)-1- 


/A(t) 


A(r)  ~ 

A(‘) = 


and 


F  = 


!rl  tl  Brl  S?  j 
!j2tl  Br2t2  i 
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C.  EMCT1QNALSXPANS1QN 

Ws  consider  here  a  method  for  reducing  the  computational  and  storage  burden  of 
the  near-field  computations.  The  basis  for  the  method  is  a  representation  of  the  measured 
data  field  by  a  suitable  expansion  function  which  is  truncated  to  limit  the  number  of  field 
function  derivatives  (to  yield  a  relatively  "smooth"  estimate  of  the  scattering  matrix)  and 
which  minimizes  some  specified  error  function  of  the  estimated  scattering  matrix.  The 
equations  are  written  in  terms  of  the  previously  determined  values,  (k,I),  of  the 
scattering  matrix  and  the  next  value,  T1:[(£-|-Ak,l+Al),  of  the  scattering  matrix  to  be 
determined.  As  an  example,  we  will  consider  an  estimate  for  the  scattering  matrix  which 
minimizes  the  norm,  N(d).  of  the  equation  defining  the  scattering  matrix,  i.e. 

N(d)  =  |  Tn(k,i)  I4rV(2r)6  jnorm. 

The  specific  matrix  normalization  which  might  be  used  need  not  be  specified  at  this  time 
for  the  purposes  of  this  example.  The  minimization  is  subject  to  the  constraint  that  in  the 
region  of  each  set  of  estimated  values,  the  scattering  matrix  has  no  derivatives  above  the 
second  order.  Therefore 

T^S+Ak.i+Ai)  =^0  (Vj-Ak  5/$k  +  y  J-Al  a/3l](n>  Tn(k,I) 
where  V  g  and  V  |  are  the  gradients  of  the  function  T,  with  respect  to  k  and  I  and 
[  ](n)  represents  the  gradient  and  partial  differential  operators  in  the  bracketed  quantity 
being  raised  to  the  nth  power  before  operating  on  the  process  functions.  The  minimization 
and  constraint  conditions  are  combined  using  a  LaGrange  multiplier  as  shown  below: 

Q  »  N(d)  +  A  |(Tn(S+Ai,I+Al)  -  1*^  [V  j-a£  8/(&  +  V  j-Al  d/dl}^ 

The  equation  for  Q  is  now  minimized  with  respect  to  A  and  also  Tj,  (£,I)  for  given  values 
of  Ak  and  Al  in  order  to  compute  Tjj(k+Ak,I+Al).  The  intent  of  the  operations  is  to 
solve  for  consecutive  values  of  f  ^(k+AS.I-i- Al)  for  spadngs  of  Ak  and  Al  which  are 
greater  than  those  which  might  be  required  if  the  process  constraints  were  not  applied.  If 
the  resultant  spadngs  of  Ak  and  Al  are  greater  than  those  that  would  be  required  based  on 
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the  effective  aperture  of  the  scattering  object,  then  we  will  have  reduced  the  density  of  data 
sampling  for  the  near-field  measurements. 

The  equations  must  next  be  written  in  terms  of  the  measurement  parameters 
(TK(k))-1,  and  (2r)6  and  the  data  variable  lf(S,I). 

The  expansion  terms  in  the  constraint  equation  are  next  written  as 

4*  d'& + 7  rAi  Tn(J>i) = + v  rAi 

a/alj(n)  (iK’(J)r1  .H(5  ,i)  •  (iRi)1)-1/^*)6- 

The  expansion  for  n=0  yields 
The  expansion  for  n=l  yields 

[V  f-Ak  5/5 k  +  V  j-Al  d/dtf1'  (!K(k)r1-F(k,I)-(IE(I)tr1-/(2x)6  = 

5/3k  iK(k) )  •  irr(k)“J.H(k,  1)  •  (ic(I)  ‘r1  + 

^^•(vg-Ak  s/aC^k,!))-^!)1)-1  + 

K(k)_1-(7  j- a!  8/dl  H(ic,I))»(IJ:(I)tr1  + 
ncdr1.H(i,I).(-i)(iE(I)t)-I.(?  J-a!  5/a!  (ird)1))  -  (n:(I)tr1]/(2x)6. 

The  expansion  for  n=2  yields 

[7  j-Ak  5/a  +  V  j-a!  a/a!j(2)  (IK'(k)r1-F(k,i)-(IE(I)lr1-/(2x)6  = 

[V2j-Ak  8*/Sk~  +  2-7  J-Ak  7  J-a!  ^/flk dl  +  72j-a!  ^/s!2]- 

oKd)  r1  •  h(£  ,  t)  -  (iu(i )  ‘r1 .  /(2w)6. 

Rather  than  complete  the  expansion  for  the  second  order  terms,  it  is  helpful  to 
consider  the  implications  of  the  expansion  tents  relative  to  reducing  the  sampling  density 
of  measurements  while  providing  a  satisfactory  estimate  of  the  scattering  matrix.  For 
example,  we  might  examine  the  computational  and  sampling  implications  of  using  a 
measurement  process  which  assured  that  the  9*TK(i)/iPi  and  $TK(i)/Sl  terms  would  be 
zero. 

The  final  steps  are  to  perform  the  minimization  with  respect  to  A  and  the 
components  of  tTTC(£)}~i -S(S,I)  *(TZ7(!)t)— 1  in  the  expansion  equations  and  solve  for 
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T11(k+Ai,I+Al).  The  estimated  value  for  T^k+Alc.I-i-Al)  is  then  compared  with 
the  measured  value  of  Tjj(J+Ajc,I-wii).  If  the  difference  between  the  two  values 
exceeds  a  specified  tolerance  value,  we  have  several  alternatives  depending  on  the  region 
over  which  the  function  minimization  was  performed.  If  the  rep  on  is  local,  the  required 
spadngs  of  the  data  samples  can  be  determined  adaptively  and  changed  as  required  at  each 
measurement  location.  If  the  minimization  region  is  global  then  the  local  error  estimates 
can  be  used  as  a  basis  for  changing  the  constraint  conditions  of  the  algorithm. 


CHAPTER  VII 

DEFINITION  OF  MEASUREMENT  LIMITATIONS 


A.  MIEOngCITOH 

The  purpose  of  this  task  was  to  investigate  the  limitations,  both  theoretical  and 
experimental,  of  the  near-fie!d  RCS  measurement  method.  We  concentrated  on  two 
topics:  (1)  data  acquisition  limitations  and  (2)  measurement  of  wide-angle  bistatic  RCS. 
Investigation  into  the  data  acquisition  limitations  considered  the  design  of  the  receiver,  the 
antenna  configuration,  and  the  storage  requirements.  The  investigation  into  wide-angle 
bistatic  RCS  concentrated  on  the  difficulties  of  performing  such  a  measurement 

Results  of  this  investigation  are  that,  depending  upon  the  frequency  of  interest, 
near— field  RCS  is  practical  with  current  technology  for  objects  ranging  in  size  from  small 
missiles  to  B-l  Bombers  for  monostatic  RCS  and  small-angle  bistatic  RCS.  For 
wide-angle  bistatic  RCS,  we  concluded  that  new  technology  and/or  modified  near-field 
theory  is  needed.  Additionally,  application  of  near-field  RCS  for  large  objects  at 
frequencies  above  1  gigahertz  is  dependent  upon  the  development  of  large  array  antennas 
appropriate  for  near-field. 

B.  DATA  ACQUISITION  LIMITATIONS 

Investigations  into  the  data  acquisition  limitations  covered  the  topics  of  receiver 
design,  antenna  requirements,  and  data  storage  requirements.  This  investigation 
demonstrates  that  near— field  RCS  is  possible  for  a  wide  class  of  problems  with  current 
technology  and  that,  with  the  advancement  of  new  technology,  it  has  the  potential  to  be 
applicable  for  many  classes  of  RCS  measurement  configurations.  In  this  section,  only 
monostatic  and  small-angle  bistatic  RCS  is  considered, 
i.  Receiver  Daiga 

A  potential  design  for  the  receiver  in  the  near-field  measurement  method  is 
shown  in  Figure  VII— 1.  The  receiver  shown  is  a  synchronous  detector  and  is  similar  to 
receivers  used  on  typical  radar  systems.  This  type  of  receiver  is  proposed  because  of  its 
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Figure  Vn-1.  Synchronous  Detector  Receiver  Design 
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Figure  VII-2.  Near-field  Measurement  Times  For  Fixed 

Aperture  Antennas 
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ability  to  rapidly  measure  the  real  and  Imaginary  parts  of  the  received  signal.  The  real 
part  of  the  RF  signal  is  obtained  by  multiplying  the  RF  signal  by  a  reference  signal  from 
the  local  oscillator  and  low-pass  filtering  the  result.  The  imaginary  part  is  obtained  by 
multiplying  the  RF  signal  by  the  quadrature  of  the  reference  signal  and  filtering  as  before. 
These  signals  are  then  converted  to  digital  form  and  stored. 

In  most  casrn,  data  storage  will  be  the  slowest  function  in  the  receiver.  However,  if  a 
buffer  is  devised  to  switch  between  several  storage  devices,  one  could  realize  a  system 
which  is  limited  by  the  A/D  conversion  and  the  settling  time  of  the  low-pass  filter.  In 
most  practical  cases,  the  A/D  conversion  time  will  be  the  limiting  factor  on  speed  if  one  is 
willing  to  use  a  wide-band  low-pass  filter  and  accept  the  increased  noise  that  wilt  be 
associated  with  the  wide  bandwidth.  Conversion  times  of  under  one  microsecond  for  12  bit 
A/D  conversion  are  common. 

Figure  VII-2  shows  the  measurement  time  required  to  make  near-field  RCS 

measurements  assuming  that  the  only  limitation  is  the  sampling  time  for  the  A/D 

conversion.  In  this  example  an  array  antenna  with  A/2  spacing  is  used,  the  A/D  conversion 

time  is  0.1  microsecond,  and  the  time  to  move  the  transmit  probe  is  assumed  negligible. 

This  graph  shows  that  for  practical  frequencies  and  for  targets  as  large  xs  the  B-l  bomber, 

A/D  conversion  times  are  not  a  limitation  with  near-field  measurements  given  enough 

receivers.  The  main  limitation  of  this  configuration,  however,  is  the  design  of  the  array 

antenna  and  recover  system,  and  the  cost  of  such  a  system.  Arrays  with  106  to  107  array 
2  3 

elements  and  10  to  10  receivers  will  be  required  at  higher  frequencies. 

2.  Data  Storage 

Data  storage  requirements  are  extremely  large  for  many  practical 
measurement  configurations.  Figure  VII— 3  shows  the  data  volume  generated  turning 
conventional  sampling  techniques  and  no  data  compression.  The  large  amounts  of  data 
that  must  be  stored  required  some  land  of  data  compression  to  make  near-field  RCS 
practical,  for  large  scattering  objects  at  high  frequencies.  Several  techniques  for  reducing 


the  volume  of  data  are  discussed  in  Chapter  VI.  Even  without  data  compression,  however, 
it  appears  feasible  to  measure  objects  ranging  from  a  jeep  at  5  Ghz  to  a  B— 1  Bomber  at 
100  Mhz. 

3.  Existing  Technology 

In  the  previous  section  it  was  shown  that  near-field  measurement  times  are 
reasonable  given  a  suitably-designed  array  antenna  and  receiver  system.  This  technology 
has  not  been  demonstrated  yet,  and  therefore  near-field  measurement  times  were 
investigated  assuming  implementations  using  existing  technology. 

From  Figure  VII-2,  note  that  measurement  times  of  weeks  are  needed  for 
large  scattering  objects  in  the  gigahertz  frequency  bands  given  an  optimally  designed 
antenna  system.  The  only  measurement  speed  limitation  of  that  system  was  the  A/D 
conversion  time.  Therefore,  one  can  conclude  that  near-field  RCS  measurements  of  large 
objects  will  require  multi-antenna  arrays  to  achieve  the  necessary  measurement  speed. 

Measurement  times  for  near-field  RCS  wee  developed  as  a  function  of  the 
number  of  probe  antennas  (or  array  elements),  the  frequency,  and  the  aperture  scan  area. 
Current  technology  generally  uses  a  waveguide  probe  antenna.  Although  this  type  of  probe 
has  limited  application  here  due  to  its  high  scattering  cross  section,  it  may  be  usable  as 
long  as  the  total  number  of  probes  is  relatively  small  (i.e.,  the  total  cross  section  ot  the 
probes  must  remain  small).  Figure  VII-4  through  VTI-8  show  the  resulting  measurement 
times  at  different  frequencies  for  various  size  apertures.  Tabk  VH-1  shows  the  practical 
limitation  on  target  size  as  a  function  of  frequency  using  these  probe  scanning  techniques. 
As  shown  in  the  table,  many  practical  systems  can  be  measured  with  current  probe 
scanning  techniques  and  an  expensive  array  antenna  system  need  not  be  developed  for 
these  cases. 

C.  MP.B-AttCLafflCTAHC.RCS 


The  near-field  theory  developed  by  Dinallo  requires  that  the  probe-to-probe 
coupling  of  the  transmitter  and  receiver  be  negligible.  This  can  be  accomplished  either  by 


Time  (day*)  Measurement  Ti 


Table  VII-1. 


Practical  Target  Size  Limitations  Versus 
Frequency  Using  Current  Probe  Scanning  Techniques. 


Frequency 

Target 

100  MHz 

B-l  bomber 

500  MHz 

F-15  fighter 

1  GHz 

M60tank 

5  GHz 

Jeep 

10  GHz 

<  2m  diameter 
target 
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using  probes  with  minimal  mutual  coupling  or  measurement  techniques  such  as  time 
gating,  FM  modulation  of  the  background  (range  filtering),  k-space  filtering  (angular 
filtering),  and  so  forth.  Each  of  these  experimental  techniques  has  its  limitations.  The 
problem  that  arises  is  that  there  will  always  be  potential  geometric  orientations  for  which 
angular,  range,  or  time  filtering  alone  cannot  distinguish  between  the  received  signals  from 
the  scatterer  and  the  transmit  probe  antenna.  This  problem  might  be  solved  by  using  two 
or  more  of  these  techniques  is  tandem.  Additional  development  of  the  theory  of  near-field 
RCS  measurement  may  also  be  fruitful.  Wide-angle  bistatic  RCS  is  conceptually  feasible, 
but  has  yet  to  be  demonstrated.  Further  investigation  into  using  multiple  filtering 
techniques  and/or  additional  theoretical  development  need  to  be  performed. 

D.  CONCLUSIONS 

This  brief  investigation  has  shown  that  near-field  monostatic  and  small-angle 
bistatic  RCS  measurements  are  practical  for  some  interesting  classes  of  problems  using 
current  technology.  It  also  shows  the  need  to  develop  both  large  array  antenna  systems 
and  inexpensive  multi-channel  receivers.  The  need  for  array  antennas,  and  thus 
multi-channel  receiving,  increases  with  increasing  target  size  and  frequency.  Data  storage 
problems  also  increase  with  increasing  target  size  and  frequency. 

The  data  quantity  reduction  techniques  discussed  in  Chapter  V  have  the  potential 
of  relaxing  the  constraints  for  both  the  data  storage  and  the  complexity  of  the  array 
antenna  by  requiring  fewer  data  samples,  and  therefore  requiring  fewer  array  elements. 


CHAPTER  VIH 
CONCLUSIONS 


This  final  report  presents  the  technical  work  performed  by  BDM  for  the  Roms  ajc 
Development  Center  to  perform  laboratory  measurements  and  analysis  to  validate  a  planar 
near-field  RCS  measurement  technique. 

The  major  issues  of  the  technical  effort  are  discussed  in  the  pra’.  ious  chapters  of  this 
report.  The  probe  corrected  near-field  theory  is  well-developed  and  is  directly  applicable 
to  the  planar  measurement  configuration.  Software  to  implement  the  theory  has  been 
developed.  Near-fiel<i  measurements  were  obtained  with  a  low-cost  computer- -oeutroUc-d 
scanning  table  designed  and  constructed  in  the  BDM  Laboratory.  Agreement  between 
measurements  and  predictions  was  excellent.  We  have  identified  several  approaches  to 
reducing  the  quantity  of  data  required  for  near-field  RCS  measurements,  and  have 
identified  some  of  the  limitations  of  the  method. 

Although  we  do  not  see  any  fundamental  obstacles  to  implementing  this  RCS 
measurement  technique  on  a  larger  scale,  our  recent  work  nas  given  us  insight  into  the  typfr 
of  development  needed  to  bring  this  concept  to  fruition.  While  building  on  the  mature 
near-field  antenna  meitsurement  technology,  planar  near-field  RCS  measurement  itself  is 
obviously  not  a  matu-c  technology.  BDM  is  not  at  a  point  to  design  turn-key  systems 
implementing  this  technology,  although  we  have  been  leaders  in  developing  and  validating 
the  mathematics  and  experimental  protocols  for  these  measurements. 

We  do  feel,  however,  that  this  Is  a  very  promising  avenue  for  near-field 
RCS  measurement.  It  shows  promise  of  addressing  several  important  issues  of  RCS 
measurement: 

(1)  measuring  RCS  of  low-observable  objects, 

(2)  measuring  monostatic  and  bistatic  RCS, 

(3)  providing  a  secure  facility  for  RCS  measurement, 


i 


(d)  increased  accuracy,  and 

(5)  increased  repeatability. 

In  the  following  sections  we  discuss  some  technology  areas  that  need  further 
development.  Specific  areas  of  planar  near-field  RCS  technology  are  discussed  in  the  terms 
of  topic  problem  areas  and  tentative  solutions.  Environmental  factors  are  also  noted  where 
appropriate. 

B.  REAHSXIC  .RCS  MEASUREMENT  REQUIREMENTS 

There  is  currently  a  credibility  gap  that  must  be  surmounted.  Although  the 
ineoretical  basis  for  -ronverting  near-field  RC>  to  far-field  RCS  is  firm,  there  has  been  no 
large-scale  demonstration  that  near-field  RCS  data  will  yield  valid,  believable  far-field 
RCS  for  realistic  military  targets.-  It  13  necessary  to  establish  analytical  bases  for  relating 
conventional  far-field  RCS  measurements  and  near-field  RCS  measurements,  identify 
limitations  of  both  far-field  and  near-field  RCS  measurements,  and  investigate  the  extent 
to  which  the  statistical  characteristics  of  fas-fidd  and  near-fiald  data  are  complementary. 

There  will  be-  oonw  difficulties  encountered  in  measuring  low-observabie  RCS,  since, 
as  in  conventional  fer-Seld  ranges,  the  RCS  may  be  comparable  to  the  background 
scattering  levels  in  the  system.  It  is  necessary  to  identify  mechanisms  for  dealing  optimally 
with  low  signal/noise  measurement  environments. 

The  actual  requirements  for  RCS  data  quantity  and  precision  are  not  clear.  While 
very  accurate  RCS  from  many  angles  is  desirable,  what  is  actually  acceptable?  There  is  a 
need  for  realistic  RCS  requirements  assessments  based  upon  system  specifications  and 
mission  requirements. 

The  analysis  of  near-field  RCS  data  has  moderate  computational  complexity. 
Studies  are  needed  to  identify  trade-off  elements,  precision,  and  computational  complexity 
and  apply  a  systems  view  of  the  problem  to  determine  computational  trade-offs. 
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C.  THEORETICAL  BASIS 

Tu2  theoretical  basis  for  near-field  RCS  measurements  needs  expanding.  Among 
the  effects  that  need  additional  investigation  are  multiple  interactions  between  probes  and 
objects,  ways  of  reducing  the  computational  burden,  effects  of  undersampling  the  data 
field,  maximum  entropy  constraints  to  reduce  what  would  be  redundant  information  in  a 
noise-ffae  environment,  algorithm  modifications  to  compute  parameters  of  estimated 
statistical  distributions  directly  for  RCS  or  antenna  patterns,  and  more  effective  probe 
antenna  designs  and  patterns. 

D.  CALIBRATION 

Existing  RF  technology  is  generally  adequate  to  support  this  RCS  measurement 
technique.  Care  must  be  taken  in  stable  design  of  signal  sources,  feeds,  and  receivers.  The 
mechanical  and  electrical  design  of  the  antennas  and  feeds  will  require  ingenuity. 
Calibration  of  such  a  near-field  RCS  measurement  facility  would  probably  require  built-in 
calibration  loops  in  the  RF  paths  to  ailow  adaptive  self-test  and  recalibration  during  data 
runs.  Implications  cf  such  self-test  loops,  and  the  possibility  of  using  created  white  (or 
other  suitably  defined)  spatial  and  frequency  noise  sources  and  out-of-band  techniques  to 
check  calibration  in  real  time,  should  be  investigated. 

E.  MEASURING  THE  NECESSARY  DATA 

The  two  key  problems  for  this  near-field  measurement  approach  are  the  amount  of 
data  that  must  be  measured,  and  the  time  that  will  be  required  to  make  the  measurements. 
There  is  no  fundamental  obstacle  to  recording  the  necessary  data,  but  it  calls  for  clever 
measurement  system  design  and  data  processing  algorithms. 

The  measured  near-field  data  will  be  voluminous.  Ways  are  needed  to  reduce  the 
amount  cf  data  required,  and  to  efficiently  store  the  data  once  they  are  measured.  Possible 
paths  to  pursue  are  to  achieve  data  density  reductions  by  alternative  computation  and 
analysis  algorithms;  data  coding  algorithms  to  reduce  data  volume  for  -tcquisition,  storage, 
and  processing;  trade-offs  between  data  precision,  storage  requirements,  and  analysis 
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accuracy;  new  algorithms  (eigecstructure,  analytic  and  series  functions,  and  orthogonal 
expansions)  for  reducing  sampling  density;  determining  what  can  be  done  numerically  to 
utilize  magnitude- -only  or  phase-only  measurements;  use  of  error-free  versus 
non-error-fret'!  data  compression;  use  of  densely-packed  antenna  arrays;  and  investigation 
of  various  hi  go-density  digital  data  storage  media  including  optical  disk  and  holographic 
storage. 

Other  avenues  of  investigation  include  fully-automated  measurements  systems, 
computer-controlled  data  collection  and  probe  and  target  positioning  with  process  quality 
feedback,  time-division  and  frequency-division  multiplex  antenna  and  feed  systems, 
fast-response  servo  control  systems  for  probe  antenna  positioning,  existing  laser  metrology 
systems  for  closed-loop  control  of  antenna  position,  multiple  transmit  and  receive  channels 
with  broad-band  intermediate,  frequency  strips  to  allow  high-speed  I-Q  measurements, 
parallel  computationaliy-sperific  dedicated  processors,  and  combined  data  measurement 
and  analog  computation. 

The  various  RF  coupling  paths  between  measurements  system  and  object  must  be 
analysed  as  part  of  the  error  analysis.  Signal-to-noise  limitations  must  be  evaluated, 
particularly  in  light  of  low-observable  targets,  broad  IF  bandwidths  in  receivers,  and 
system  noise  figure.  Other  signal-to-noise  issues  include  using  multiple  coherent  data 
measurements,  background  scatter  reduction,  coherent  or  incoherent  background 
subtraction  and  suppression,  absorption  and  signature  techniques,  effects  of  multiple 
reflections,  explicit  time  gating  and  implicit  gating  using  synthetic  spectra  and  data 
signatures,  environmental  changes  and  deviations  from  Ideal,  and  factor  analysis. 

Part  of  that  investigation  would  be  to  evaluate  the  level  of  specular  and  diffuse 
reflections  that  may  be  present  from  both  object  and  background,  and  possibly  create 
modified  algorithms  tailored  respectively  to  emphasise  or  reject  specular  and  diffuse 
reflections  via  a  feed-back  loop,  using  results  to  isolate  undesired  components  of  the 
background. 


Obtaining  the  necessary  metrology  accuracy  is  not  a  trivial  problem,  but  is  well 
within  the  state  of  the  art.  Tools  to  be  used  here  include  laser  interferometry,  closed  loop 
control,  real-time  RF  frequency  and  phase  error  measurement  and  adaptive  correction,  and 
phase  locking  techniques.  Timing  errors  can  be  minimised  by  referring  all  measurements 
to  a  master  clock  to  correct  the  time  references  of  data  as  required. 

It  is  necessary  to  evaluate  reproducibility  of  measurement  conditions. 
Reproducibility  is  impacted  by  environmental  effects,  equipmait  drift,  structural  stability, 
target  stability,  and  aging,  and  might  be  measured  in  terms  of  pattern  distributions  and 
statistical  confidence  intervals. 

Effects  of  non-ideal  measurement  components  (e.g.,  amplifiers,  mixers)  must  be 
evaluated.  One  approach  is  to  exercise  and  modify  algorithms  as  required  to  examine  data 
behavior  given  non-ideal  components,  including  error  analyses,  simulation,  and 
experiments  to  evaluate  effects  of  component  variability. 

Possible  methods  of  alleviating  undesired  effects  of  probe  characteristics  include 
using  broad-band,  low  mutual  coupling,  variable  polarization,  active  nonlinear  elements. 
Adequate  polarization  and  geometry  control  may  allow  simplifications  in  the  data 
processing  algorithms;  this  will  require  detailed  antenna  analysis  and  design  and  feed 
systems  using  PIN  diode  switching  and  compensating  elements. 

Environmental  noise  components  might  be  isolated  and  compensated  for  by 
measuring  and  assigning  data  signatures  to  various  environmental  and  measurement 
elements  and  demodulating  the  undesired  signatures  in  real  time.  Environmental 
interference  might  also  be  reduced  by  using  out-of-band  data  from  known  targets  to 
isolate  environmental  interference  on  a  real-time  basis. 

Approaches  to  reducing  data  computational  burden  indude  using  canonic 
simplifications  based  on  alternative  measurement  and  analysis  algorithms,  and 
simultaneous  measurement,  computation,  and  analysis  using  buffered  computations  and 
direct  transform  algorithms. 
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Data  collection  monitoring  and  quality  control  call  for  automated  operations  and 
on-line  factor  analysis  to  determine  data  error  sensitivity.  The  need  for  Identifying 
unusual  measurement  problems  during  data  acquisition  suggests  that  one  should  create  an 
evolving  data  pattern  generated  incrementally  as  measurement  data  are  obtained,  so  that 
algorithms  can  detect  divergence  from  the  emerging  pattern. 

F.  DATA  ANALYSIS 

The  relations  between  input  and  output  data  quantity  and  quality,  and  numerical 
effects  in  the  data  analysis  algorithms,  need  additional  investigation.  Possible  approaches 
are  to  employ  a  64-bit  (or  longer)  word-length  computer  to  simulate  a  variable  word 
length  computer  and  thus  determine  the  numerical  sensitivity  of  the  algorithms  and 
optimize  the  algorithms  for  minimum  sensitivity;  investigate  dynamic  range  limitations 
due  to  numerical  considerations  of  using  large  Fast  Fourier  Transforms  (FFTs);  and 
possible  use  of  phase-only  data.  Simulations  are  needed  to  help  determine  error  relations 
between  data  inputs  and  outputs.  Data  analysis  time  and  resources  could  be  reduced  by 
using  dedicated  hardware  FFTs  and  processors. 

Additional  paths  to  reducing  computational  complexity  include  using  the  physical 
situation's  constraints  (both  target  and  environment),  using  mathematical  structure 
constraints,  adding  other  constraints  on  information  content,  matching  the  information 
content  of  the  channel  between  measurements  and  canonic  results,  factor  analysis,  and 
various  operations  research  techniques  (including  linear  programming). 

Other  possible  data  and  algorithm  modifications  include  using  blank  data  cells 
(under  sampling),  truncation  (blanking),  zero  interpolation,  blanking  interpolation,  various 
simplifications  and  alternatives  to  the  existing  near-field  theory,  solving  dual/analog 
problems,  Identifying  relations  between  simplification  and  useful  outputs,  computing 
pattern  statistics  rather  than  RCS  patterns,  and  using  object  and  algorithm  functional 
symmetry. 
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G.  OPERATIONAL  CONSIDERATIONS 

We  emphasize  again  that  it  is  not  dear  what  RCS  data  are  really  needed  for 
existing  aircraft  and  missiles.  It  would  be  a  noteworthy  capability  simply  to  be  able  to 
measure  monostatic  RCS  at,  say,  one  degree  intervals  for  angles  within  *30  degrees  of  the 
horizon,  but  it  is  not  dear  whether  mission  effectiveness  would  be  truly  enhanced  by  such 
detailed  knowledge.  Perhaps  it  would  be  more  useful  to  be  able  to  state  that  the  maximum 
RCS  over  a  specified  solid  angle  is  less  than  s  with  probability  y.  The  range  of  angles  and 
the  angular  resolution  of  the  data  have  major  technical  (measurement  time  and  data 
volume)  and  economic  (time  and  cost)  impact.  As  aids  to  baselining  this  measurement 
approach,  time  and  data  volume  could  be  estimated  for  specific  measurement  scenarios, 
given  the  angular  range  and  resolution  needed. 

The  absolute  RCS  accuracy  required  for  mission  evaluation  has  a  major  impact 
upon  how  much  data  are  needed,  how  long  it  takes  to  measure  the  data,  and  what  kind  of 
algorithms  are  used  to  process  and  display  the  data.  Acceptable  and  desirable 
measurement  accuracy  constraints  for  mission  scenarios  must  be  established. 

The  mission  scenario  elements  together  impose  stringent  requirements  upon  the 
measurement  system.  As  part  of  identifying  true  measurement  needs,  a  system  level 
evaluation  of  mission  and  operations  use  of  RCS  data  is  needed,  to  determine  feasible  costs 
(time,  capital,  and  staff),  and  to  perform  trade-off  analyses  and  simulations, 
i 
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APPENDIX  A 
COORDINATE  SYSTEMS 


A.  COORDINATE  SYSTEMS 


There  are  two  coordinate  systems  (figure  A — 1)  to  deal  with  in  performing 
near-field  measurements:  the  scan  table  coordinate  system  (S')  and  the  antenna 
coordinate  system  (S).  In  the  context  of  bistatic  RCS  measurements,  the  antenna 
coordinate  system  is  referred  to  as  the  test  object  coordinate  system. 

The  scan  table  coordinate  system  is  the  "laboratory"  coordinate  frame  in  which 
measurements  are  taken,  while  the  antenna  coordinate  system  is  the  natural  coordinate 
frame  for  the  antenna  or  object  under  test.  Both  are  right-handed  coordinate  systems. 

The  scan  table  coordinate  system  is  defined  with  horizontal  x'—y'  plane  and 
z '-axis  upwards.  The  x'-axis  corresponds  to  movement  of  a  platform  only,  while  the 
y'-auds  corresponds  to  movement  of  a  carriage  along  with  its  platform. 

The  origin  of  the  scan  table  coordinate  system  is  at  the  center  of  the  aperture  of  the 
TX  probe  antenna  when  the  TX  platform  is  at  its  initial  ("home")  position.  The 
horizontal  position  of  the  origin  coincides  with  the  TX  probe  antenna  mounting  hole  in  the 
TX  platform,  while  the  height  of  the  origin  above  laboratory  floor  level  depends  upon  the 
height  of  the  TX  probe  antenna  that  is  installed. 

The  origin  of  the  antenna  coordinate  system  S  is  placed  at  some  convenient 
reference  point  on  the  AUT.  If  r'  is  the  representation  in  S'  of  a  vector  from  the  origin 
of  S'  to  the  origin  of  S,  then  a  vector  r'  expressed  in  S'  transforms  to  S  as 
?  =  ¥-(r'-r^), 

where 


ro5  roxex  +  roy®y  +  roz®z’ 
rl  0  (T 
0-1  0 
.0  0-1. 

rsrx5x+ry6y  +  rzVaad 
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If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  as  (r ,8,4>),  then 
its  negative  (—  A)  is  written  as  (r,r—  8,x  +  41)- 

If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  in  S  as  (r,M)> 
then  it  is  expressed  in  S'  as  A'  *  (r',$',^'),  where 
r'  =r 
S'  =  1-8 
—~4>- 

If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  in  S  S3  (r ,8y<j>), 
then  its  negative  (—A)  expressed  in  S'  is  —A  *  (r'.tf',^'),  where 
r'  =  r 
8’  =  8 
—  -{*  +  <!>)• 

The  above  relations  can  be  verified  by  inspection  of  the  coordinate  systems 
involved.  We  require  that 
r  >  0 
0  <8<x 
—x<4><x. 

S'  is  used  to  describe  the  probe  antenna  TX  and  RX  characteristics  since  these 
functions  do  not  depend  on  a  probe's  location  in  S'. 

ITJJIMINATIFG  PLANE  WAVE 

The  target  is  illuminated  by  a  plane  wave  with  polarization  (£j*k£=0)  that  is 
effectively  emitted  from  coordinates  (xj>y^{)  in  S'  toward  the  target  centered  at  known 
r<YM7«n»tP«  (x^,y',z£)  in  S'.  Define  a  unit  vector  rj  in  S'  in  the  direction  (x^-xj,  yjj-yp 
s^-zj).  The  vector  rj  is  known,  since  rj  is  a  unit  wave-vector  from  the  far-field 
illumination  direction  (if  TX  probe  is  scanning),  or  (xj  .y^)  is  the  TX  probe  position  and 
r'  can  be  calculated'  (if  the  TX  probe  is  stationary). 
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If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  as  (r ,8,4>),  then 
its  negative  (—  A)  is  written  as  (r,r—  8,x  +  41)- 

If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  in  S  as  (r,M)> 
then  it  is  expressed  in  S'  as  A'  *  (r',$',^'),  where 
r'  =r 
S'  =  1-8 
—~4>- 

If  the  unit  direction  vector  A  is  expressed  in  spherical  coordinates  in  S  S3  (r ,8y<j>), 
then  its  negative  (—A)  expressed  in  S'  is  —A  *  (r'.tf',^'),  where 
r'  =  r 
8’  =  8 
—  -{*  +  <!>)• 

The  above  relations  can  be  verified  by  inspection  of  the  coordinate  systems 
involved.  We  require  that 
r  >  0 
0  <8<x 
—x<4><x. 

S'  is  used  to  describe  the  probe  antenna  TX  and  RX  characteristics  since  these 
functions  do  not  depend  on  a  probe's  location  in  S'. 

ITJJIMINATIFG  PLANE  WAVE 

The  target  is  illuminated  by  a  plane  wave  with  polarization  (£j*k£=0)  that  is 
effectively  emitted  from  coordinates  (xj>y^{)  in  S'  toward  the  target  centered  at  known 
r<YM7«n»tP«  (x^,y',z£)  in  S'.  Define  a  unit  vector  rj  in  S'  in  the  direction  (x^-xj,  yjj-yp 
s^-zj).  The  vector  rj  is  known,  since  rj  is  a  unit  wave-vector  from  the  far-field 
illumination  direction  (if  TX  probe  is  scanning),  or  (xj  .y^)  is  the  TX  probe  position  and 
r'  can  be  calculated'  (if  the  TX  probe  is  stationary). 
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The  propagation  vector  &£  is  expressed  as 
kil  =  rilk 

kyl  =  rylk 

kzl  “  rzlk 
or 

rl=k 

«i=cos  \t'z1) 

<t>{  =  tan-1  (r^pr^) 

where 

ks2x/A. 

We  also  know,  or  will  specify  arbitrarily,  the  polarization  of  the  illuminating  plane 
wave  If  the  TX  probe  is  scanning,  we  select  the  polarization(s)  Eg  or  E^  for  a  particular 
simulation.  If  the  TX  probe  is  a  stationary  standard  gain  horn  pointed  at  the  object,  we 
will  orient  the  horn  so  that  the  object  is  on-axis  and  the  E-plane  of  the  horn  is  either 
vertical  or  horizontal  in  S',  so  that  the  transmitted  polarization  is  either  purely  E'g  or 
purely 

If  the  polarization  is  known  in  the  basis  set 

then  the  spherical  components  of  the  polarization  are 

E'j  =  E^jSinfljCOS^j  +  E^sin^sin^,  +  E^coeflj 
~  E^icos^cos^i  +  E^cos^sin^,  -  E'jrintfj 
~  -E^sin^i  +  E'jcos^. 

If  the  polarization  is  known  in  the  basis  set 
Ej  e  (E'^jj^j), 

then  the  Cartesian  components  of  polarization  are 

E^j  =  E'jSinSjCoe^j  +  Ejjcoefljcos^j  -  E^sin^j 
Eyl  “  E'fSin^sin^i  +  E^cos^sin^j  +  E^cos^j 
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Ezl  =  EnC03dl_Etfl8inei> 

and  E'  must  be  selected  to  ensure  that  Thus  can  be  expressed  in  either 

Cartesian  or  spherical  coordinates  using  the  above  transformation  equations. 

The  propagation  vector  is  expressed  in  S  via  the  transformations 

kxlas+kil 

kyl  ~  "^yl 

kzl  =  ^l 
or 


ri=r'i 

0i  =  — 0V 


The  polarization  vector  is  expressed  in  S  via  the  transformations 


Exl  =  +Eil 
Eyl  =  ~Eyl 
Ezl^l 


or 


E«i = ~Eh 
E0i =  ~Ek- 

Therefore  we  can  express  the  incident  wave-vector  (k^)  and  polarization  (£j)  in  either 
Cartesian  or  spherical  components  in  S,  which  is  the  coordinate  system  of  the  target. 

D.  SCATTERED  PLANS  WAVE 

The  target  scatters  the  illuminating  plane  wave  into  plane  waves  in  various 
directions.  In  order  to  evaluate  the  scattering,  we  select  a  particular  scattering  direction 
based  upon  where  the  observer  is  in  the  far-field. 

Consider  a  plane  wave  that  is  scattered  toward  the  known  coordinates  in 

S'.  Define  a  unit  vector  r£  in  S'  in  the  direction  (x^-Xg^-ygiZ^-ZQ).  The  propagation 


vector  kj  is  expressed  as 

•  kJa!“ri2k 
ky2s=rf2k 

kz2  =  tz£ 
or 


where  k  s  2r/A. 

If  scattering  in  a  known  direction  is  desired  rather  than  scattering  toward  a  known 
coordinate,  then  (expressed  in  S')  is  the  scattering  direction  and  (x^y^.Zj)  need  001  be 
known. 

The  propagation  vector  k^  is  expressed  in  S  via  the  transformations 

kx2  =  +kx2 
ky2  =  -ky2 


At  this  point  we  know  the  incident  k-vector  (kj),  the  incident  polarization  (E^), 
and  the  scattering  k-vector  (kj),  all  expressed  in  S,  so  we  can  compute  the  polarization 
vector  (fij)  of  the  far-Geld  scattered  plane  wave  (Figure  A— 2)  by  means  of  an  analytic  or 
numerical  model  of  the  scattering  object.  The  usual  range-dependent  phase  factor  can  be 
ingored. 
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Scattering  Object 


I!lum>n*t!on 

from 

Far-FIeld 

(  8 j  ) 


Observer 

in 

Far-Field 

(82**2) 


Figure  A-Z  Scattering  Geometry 


The  scattering  analysis  of  the  sphere  requires  that  the  object  be  .  in  a  preferred 
coordinate  system  S"  in  which  it  is  illuminated  by  an  incident  plane  wave  having 
propagation  vector  |kjez  and  linear  polarization  jS|ex.  Since  kj  and  are  not  so 
constrained,  we  must  construct  a  coordinate  rotation  matrix  £  that  defines  a  Dew 
coordinate  system  S'1  in  which 
A-itj  =  IkjjejJsnd 

where  we  must  have 

V®2=0- 

The  coordinate  rotation  A  transforms  an  arbitrary  linearly  polarized  incident  plane 
wave  into  a  plane  wave  propagating  in  the  e“  direction  with  e£  polarization.  The 
scattering  direction  is  also  transformed  as 

At  this  paint  kp  and  kj  are  known  numerically,  so  the  scattered  polarization  £lj 
can  be  calculated  in  S"  using  the  appropriate  analysis  (e.g.,  Stratton  for  the  dielectric 
sphere)  and  then  transformed  back  to  S  as 

The  matrix  A  is  computed  in  three  steps.  Let  Ax  be  a  matrix  that  rotates  a  vector 
k  (for  which  we  will  actually  use  kj)  about  §x  by  an  angle 

—1  ky 

O’.  =  —  tan  -rf—  . 

z 

This  rotation  reduces  ky  to  zero,  so  the  rotated  k  (which  we  will  call  k^)  lies  in  the  exez 

plane.  The  matrix  is 

1  0  0 
Ax=  0  cos^j)  sinfaj) 

0  -si^Oj)  cos(a1) 

so 
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Ed), 


■A** 


and  applying  the  same  rotation  to  fe  gives 

fed)  =  ^.fe. 

Next,  let  Ay  be  a  matrix  that  rotates  a  rector  about  ey  by  an  angle 

kd>i 


■  -tan 


-1 


x 

vd) 

••  “z 


This  rotation  reduces  k£d  to  aero,  so  the  rotated  £d)  (which  we  will  call  fed))  lies  on  the 


e  axis.  The  matrix  for  this  rotation  is 
z 


cos(o^)  0  sin^) 
0  1  0 
-sin(o^)  0  cos(o^) 


so 


fed)  *  Ax-£d) 
fed)  =  ax-^(1). 

At  this  point,  fed)  ^  on  the  ez  axis,  as  desired,  and  fed)  is  in  the  exey  plane. 


Finally,  let  A_  be  a  matrix  that  rotates  a  vector  about  e  by  an  angle 


“3 


=  tan 


-1 


E^2)] 

-* 


This  rotation  reduces  to  zero,  so  the  rotated  fed)  (which  we  will  call  fed))  will  lie  on 
the  §x  axis.  The  matrix  is 


008(03)  sin(og)  O' 

Aj  *  |  -sin(o^)  cos(c^)  0  . 

0  0  1, 

The  desired  coordinate  rotation  matrix  is 


A  ”  ^2 ' 4y '  — x' 

The  scattering  analysis  of  the  disk  requires  that  the  object  be  in  a  preferred 
coordinate  system  SK  in  which  it  is  illuminated  by  an  incident  plane  wave  from  above  the 
e  H-e  "  plane  having  propagation  vector  fe  =  (-ein0"e  -  cos0"ez)|k|  and  linear 

*  y 


polarization  E^ejj  +  E^e^,  where  ejj  Is  in  the  plane  k-e“  and  e^aSy.  Since  ^  and  Ej 
are  not  so  constrained,  we  most  construct  a  coordinate  rotation  matrix  A  that  defines  a 
new  coordinate  system  S"  in  which 

A-kj  =  l^j  (-sia5"§x  - cos^'ljj)  and 

aV>W%A 

where  we  must  still  have 


VV°- 

The  scattering  direction  is  also  transformed  as  k^  =  A^- 

The  matrix  A  is  computed  in  two  steps.  Let  Ax  be  a  matrix  that  rotates  a  vector  k 


(for  which  we  will  actually  use  kj)  about  ex  by  an  angle 


0  if  kz  <  0 

r  if  k  >0 
2 


This  rotation  ensures  that  is  non-positive,  as  required,  so  the  rotated  k  (which  we  will 

call  k^)  lies  below  the  e  e  plane.  The  matrix  is 
*  y 

1  0  0 
Ax=  0  cosfap  sinfaj)  , 

0  -sin(aj)  cos(o^) 


so 


and  applying  the  same  rotation  to  £  gives 
£^  =  Ax-E.  • 

Then  let  A2  be  a  matrix  that  rotates  a  vector  about  e£  by  an  angle 


Thi3  rotation  reduces  to  zero  and  makes  kj^  negative.  The  matrix  is 
cos^)  sin(o^)  0 
Az=  -sin(o^)  cos(o^)  0  . 

0  0  1. 
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■;-r 
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The  desired  coordinate  rotat:  on  matrix  is  £  =  4s'4x‘  Applying  ^  to  £p 

^1  =  4^ 

=sSlx6x  +  Ely6y  +  E  lA 

In  S"  the  angles  of  incidence  and  emergence  of  the  plane  waves  are 

q  ~  coT1  (Icy kj) 

rfJ-C, 


^-^(kyy 


=  tan-1 


The  incident  polarization  parallel  and  perpendicular  to  the  plane  of  incidence  in  S"  are 
%j-  I^K-eosSJeJ  +  dnflJep 

At  this  point  kj,  Sp  and  k£  are  known  numerically,  so  the  scattered  polarization  Eg  can  be 
calculated  in  S"  as  a  linear  combination  of  E%  and  E^: 
r  BSh  1  rs„  S10UE?„1 


311  ^12 1  [El||] 

321  S22  UlJ 


In  terms  of  egegeg. 


^  =  E2ir¥x+«  +  E21J^z- 


“2  —  2|!“”p2  x  "r  °2iey  ‘ 
which  is  transformed  back  to  S  as 


*8 


M  r J  fttlO  *  u  W  1 W '  I*)  mVfViH 


If  the  predicted  fax-field  scattered  plane  wave's  polarization  vector  is 
h  s  when  expressed  in  S,  then  the  spherical  components  of  polarisation  are 

Erk  =  ExjSinfijeos^  +  E^slnfc^in^  +  E^cos^ 

E02  =  E^cos^cos^o  +  EyjCOsSjSin^j  -  E^in^ 

E*B"'*x2?to*2  •fEy2c035>2- 
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If  the  calculated  polarization  in  spherical  coordinates  is 
®2  3  (Er2>Effi>E#) 

when  expressed  in  S,  then  the  Cartesian  components  are 
Ex2  =  Ejgsinfyse^g  +  E^cos^cos^  ~ 

Ey2  =  Ej^sinffjSln^  +  E^cos^npj  +  E^cos^ 

E22  =  Er2c03^2~E52sia52- 

The  polarization  vector  is  expressed  in  S'  via  the  transformations 


Erf 

=  +Ex2 

“-Ey2 

Sz2 

=  _Ez2 

e;2 

=  +Er2 

E52 

E*2 

By  means  of  the  above  equations,  the  scattered  wave-vector  (k^)  and  polarization 
<%>  can  be  expressed  in  S'  using  either  Cartesian  or  spherical  components. 

The  fax-field  scattered  electric  field  is  given  by  fig  and  kg,  which  can  be  expressed 
in  terms  of  S  or  S'.  It  is  probably  best  to  stay  in  S,  which  i3  the  natural  coordinate  system 
for  the  scattering  object. 
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B. 


In  this  section  we  summarfce  the  relationship  between  the  angular  plane-wave 

spectrum  cf  an  B-fie!d  and  the  far— field  E-field  that  emerges  by  applying  the  method  of 

stationary  phase  to  calculation  of  fir-field  E-field. 

Consider  an  E-field  expressed  in  the  following  form 

,,  ff, ,  i& 

E{r)  =  jJ  B(k)  e  dkxdky. 

The  first  step  in  finding  the  far  field  is  to  convert  to  spherical  coordinates.  Thus, 
we  define 

P (6,4>)  s  k“sin?cos?  B(k) 

5  s  k  (sin?  cos$  ex  +  sin?  sin^  ey  +  cos?  §2) 
r  3  r  (sin®  cos£  ex  ~  sin®  jin^  ey  +  cos®  e2), 

from  which  for  kjjr  we  obtain 

ikrtcos?  cos®  +  sin?  sin®  cos(f®}] 
e  d?dtf. 

This  »  now  in  the  correct  form  for  invoking  the  method  of  stationary  phase  [3]. 

.  Since  the  method  is  somewhat  involved,  we  quote  the  result  and  reference  the  details  to 
Born  and  Wolf.  In  their  notation  (Appendix  HI,  pp.  753-754,  eq  16  and  19), 


B-3 


and  (xQ)y0)  is  the  critical  point  at  which 

flf  _  di  . 

*T5c - 3T“°- 

Thus,  is  oar  notation,  the  far  field  is 

ikr 

limit  fe{?)  -  — coe<fl)  g(£) 


r  -w 


Ca 


2£|~cos(#)  [B/tJ)e,+  B 


where 


Suppose  we  have  a  far  field 


ikr 


jimi  t  £(r)  *  a0[E^(J,^)e^  +  £^C*,$e^j 
Substituting  into  the  previous  equation  pve3 


i  a„ 


and 


2*k  ccs# 


i  a. 


B/W)=nrr¥f'^ 

y  2*k  cos#  v 

Thu *  the  acgdsr  plane-wave  spcctimn  §(k)  of  an  E-fieid  can  be  expressed  in 
terms  of  the  known  far-Celd  E-ficJd  by  the  relation 

where  kj|r. 

Later,  when  using  plane  wave  spectra  to  define  the  transmit  characteristic  cf  an 
antenna,  we  use  the  symbol  s1Q(k)  instead  of  B(k).  Hence, 
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C.  OPEN-ENDED  WAVEGUIDE 

In  this  section  we  determine  the  magnitude  of  the  TEjq  mode  excited  in  an 
open-ended  waveguide  (OEWG)  when  illuminated  by  a  plane  wave.  The  motivation  for 
this  analysis  is  clarification  of  the  field  relations  and  coefficients  of  proportionality. 

The  geometry  of  this  problem  is  illustrated  in  Figure  B— 2.  Assume  that  the 
incident  E-field  at  the  origin  is  a  plane  wave 
Etoc-E^e(W)ej+  E }nC(W)e^, 

where  E^nc(0,0)  and  E^nc(0,$  are  complex  scalars  with  an  identical  phase. 
We  will  create  an  elemental  dipole  in  the  far  field  that  could  excite  this  plane  wave. 

The  E-field  of  a  dipole  is  [4] 
ikr 

E^Pfr)  «  ■  ■'fo—  sin(tf)  IdT, 

-iij  k2 

where  we  have  defined  k  *  — ^-jj — ,  and  i?  is  the  intrinsic  impedance  of  free  space. 

We  first  choose  a  dipole  in  the  far  field  parallel  to  ^nc  to  create  E*nc  at  the  origin: 
ikr 

(Ml)* 

ikr 

where  dl  is  parallel  to  6*°°.  Note  that,  in  the  sin(0)  factor  in  £^p(r),  ff  is  the  angle 
between  dl  and  k  and  is  therefore  x/2. 


So, 

Idl=-JiFEinc, 

k  e 

which  produces  the  observed  incident  field  S*20. 

We  now  wish  to  find  the  wave  launched  into  the  OEWG  by  the  dipole  IdT.  To  do 
so,  we  invoke  reciprocity  and  calculate  the  E-field  at  the  dipole  due  to  an  excitation  of  the 
OEWG.  Suppose,  in  the  waveguide,  we  launch  a  wave  of  magnitude  E0  towards  the  end  of 
the  OEWG.  At  the  aperture,  we  have 
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figure  B-2.  Open-End  Weveguide  Coordinate 

System  for  Fir  field 


Ey  =  Eo  (1+r)  005  [tt]  »“d 
Hx=— ^-d-r)«»[¥] 

where 


Next  we  determine  &  current  distribution  that  excites  this  mode  in  the  OEWG  at  a 

point  A  from  the  opening  (Figure  B-3) 

© 


where 


V 


~T  *■  is 

1  \E~aj 


Hence,  the  above  current  produces  the  open-ended  waveguide  far  field 
ikr 

£(r)  =  A£  , 


where  for  notational  convenience  we  have  defined 


and  Eg(0,0),  Eg(tf,0),  and  are  as  defined  by  Yaghjian.  Note  that  and 
simple  unitless  functions. 

If  we  now  generalize  for  an  arbitrary  antenna,  we  find 


f(W)  -  Umit  • 

The  reciprocity  theorem  states  that 

JJJ  iMb  dv  =  jjj  6b-3a  dv. 


are  now 
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Figure  B-3.  Open-End  Waveguide  Reference  Plane 
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In  this  instance  the  two  cases  (a)  and  (b)  are: 

(a)  the  current  sheet  in  the  OEWG  excites  the  far  field,  and 

(b)  the  small  dipole  in  the  far  field  excites  the  TEjg  mode  in  the 
OEWG. 


Now,  the  mode  excited  by  the  plane  wave  in  the  OEWG  at  a  distance  inside  the 


open  end  is 

Ey  = 

where  bQ  is  unknown.  Invoking  reciprocity,  we  find 


Ey  =  boC0S(^-) 


i>-2.  ._„[¥]V 


ikr 


gl  =  AEJF[i?e*«+f5%]’ 

jb  =  idT, 


and 


Eb=bocos[¥]v 

The  left-hand  side  of  the  reciprocity  equation  above  becomes 

JJJe*-3bdv_  aev[^, + 1 $%] 

:  +  l 


while  the  right-hand  side  becomes 

f+b/2  r+a/2 


6b-JadV  = 


— b/2 


-2boEo 


t 


dx  dy 


r+b/2  r+a/2 


-b/2 


cos 

—a/2 


w 


:|  dxdy. 
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Using 


J-a/2 
the  integral  is 

JJj£b-3adV= 

Equating  the  left-  and  right-hand  sides  of  the  reciprocity  equation, 


Using  the  relation 


The  above  equation  gives  us  the  magnitude  bQ  of  the  TEjq  mode  excited  in  an 
open-ended  waveguide  when  illuminated  by  an  incident  field  Ebic.  The  receiving 
characteristic  of  the  OEWG  is  thus 


The  ratio  AE/EQ  is  calculated  elsewhere. 

Now  let  open-ended  waveguide  (OEWG)  be  illuminated  by  an  antenna  with  known 
far— field  pattern  and  excitation.  The  waveguide  feed  for  the  antenna  is  assumed  to  be 
identical  to  the  waveguide  used  as  the  OEWG.  The  geometry  for  this  analysis  is  shown  in 
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Figure  B-4.  Note  that  the  source  antenna  is  placed  in  the  unprimed  coordinate  system, 
and  the  OEWG  probe  is  in  the  primed  system.  We  define 
aQ  =  excitation  in  waveguide  of  known  antenna,  and 
bQ  =  excitation  in  OEWG. 

Both  aQ  and  bQ  are  defined  at  a  distance  from  the  origin  of  the  respective  coordinate 
systems. 

In  the  coordinate  system  0  (that  of  the  antenna  under  test),  the  far-field  E-fieid  of 
the  antenna  is 

E(f)  =  a0  [E^(^)etf  +  E^(tf,^)e^  +  0-efj  Sgj-, 
where  E and  E^(0,<6)  are  known  functions,  aQ  is  the  excitation,  and  angles  0  and  <t> 
are  related  to  the  propagation  vector  k  by 


It  has  been  shown  that  E(r)  can  be  expressed  in  terms  of  the  Fourier  transform  of 
the  components  of  £(r)  that  lie  in  the  x-y  plane.  The  Fourier  transform  of  those 
components  is  commonly  referred  to  as  the  angular  spectrum  (of  plane  waves)  of  E(r)  and 
is  denoted  here  by  B(k).  Using  the  method  of  stationary  phase,  the  angular  spectrum  of 
E(r)  can  be  related  to  the  far-field  El— field,  making  the  above  expression  for  E(r) 
equivalent  to 

£(r)  =  5(k)  e^*r  dkxdky, 

where 

S(£)“i^Nw'+E*w*)8*l' 
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Source  Antenna 


Figure  B-4.  Coordinate  Conventions  fcr  Horn  Measurements 
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or,  in  more  general  notation, 


sio(k) =  ‘ 


•limit  kre' 


j 


itk  cost)  jr  |  -vs  : 

These  plane  wayes  propagate  from  the  antenna  to  the  origin  O'  of  the  primed 

ik  d 

coordinate  system.  When  they  reach  O',  they  are  shifted  in  phase  by  a  factor  e  ,  so  the 
field  transmitted  to  O',  expressed  ja  terms  of  coordinate  system  0,  is 
ia„ 


E(0,0,d)  = - T 

2rir 


COS{D) 


ik  d 

e  2  dicxdky. 


If  now  the  probe  is  scanned  over  the  plane  z=d,  we  introduce  an  additional  phase  factor  to 
get 


E(x,y,d)  =  — ^ 
2xk‘ 


j  +  E, 

'  - mjf 


i(k  x+k  v)  ik  d 


>  e  2  dk  dk  . 

x  y 


Having  established  the  incident  field  at  the  OEWG  as  a  summation  of  plane  waves, 
the  naxt  step  is  to  calculate  the  OEWG  excitation  bQ  due  to  a  single  incident  plane  wave. 
Consider  a  plane  wave 

E(J')  =  [E/ff'^'Je^+E^tf'^Oe^k'-k'^kp-k^) 

with  propagation  vector  k',  expressed  ia  terms  jjf  coordinate  system  O',  incident  upon  an 
open-ended  waveguide.  It  was  shown  above  that 

»« - — ■*■>+ %(>■■*%(•■.*■)] 

where  I'JJ' ,tf>')  and  are  the  known  pattern  factors  of  the  OEWG: 

®OEWG^  =  AE  TcT  f|jJ(^',?,)e^]. 

In  fact,  a  plane  wave  from  the  antenna  at  0,  and  incident  upon  the  OEWG  at  O', 

in  the  unprimed  coordinated  of  the  antenna  under  test,  is 
ik.d 
aoe 


ia„e  '  i(k  x+k  v) 


2sk^cosd 
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ik  d 


30 


— 4ri 


ikd 

“oe  i(^y)r 


abk4|  1  -  |j£jJ2  2jtco«(0) 
Simplifying,  and  using  the  relations 

where  9‘—9  and  $'=*>-0,  we  get 


*'A')  + 


bo  = 


AE 

■b“ET 


abk* 


4j  l-J^f  coaW 

Summing  over  a  continuum  of  plane  waves: 


.  _AE 

b°~V 


-2ia„ 


abk 


,i.4 n  r  li 

J  IFa) 


E  ,(9,0^,*-+)  + 

coeXTj  K 
ik2d  K^jjX+kyy) 


•e  "  '  dkxdky. 

The  above  is  a  complete  and  rigorous  description  In  the  coordinate  system  of  the 
antenna  under  test.  In  general,  however,  the  measurements  will  be  made  in  a  coordinate 
system  that  is  attached  to  a  fixed  point  on  the  scanning  table,  and  is  thus  defined  by 


V 

1 

0 

o' 

x-x0' 

y" 

SS 

0 

-1 

0 

• 

y-y0 

.z" 

.0 

0 

-1 

z-d 

where  the  origin  of  the  coordinate  system  is  at  the  probe  antenna  coordinates 
0"  =  (x0,y0,d).  Hence,  after  taking  the  data,  the  first  step  of  the  analysis  is  to  convert 
from  the  doubly-primed  coordinate  system  of  the  scanning  table  to  the  unprimed 
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coordinate  3/stem  of  the  antenna  under  teat.  We  now  have 


The  transmit  characteristic  of  the  probe  antenna  is  related  to  its  far-fidd  E-Md 


sin(2)  =— i - limit  •  kre"^'^l(rl}, 

10  2*k2cae0  rn{  J 


where  k  is  parallel  to  r.  Assuming  a  reciprocal  antenna,  the  receiving  characteristic  is 
- 

roi(k; - T0 - sioHO- 

The  constants  YQ  and  ij0  are  admittances  that  may  be  nounslteed  to  anity.  If  we 
wish  to  perform  an  absolute  measurement,  it  will  be  necessary  to  specify  the  ratio  Yc/:?0 
exactly.  To  do  so,  we  return  to  the  expression  for  the  received  signal  derived  earlier 


ikzd 


have  therefore 


— 2iae  ae  ff  tjiwjA*-*)  CV 


abk*  1- 


"TT2~Eo“ 


Ky+y) 


dkxdky. 


where  bQ,  Ag,  and  EQ  are  all  in  units  of  volts/ meter.  By  equating  the  two  expressions  we 


-2*i  rAEl  M)*M 


FtM 


M4B1 

cose 


From  a  previous  section  we  have 

^iq(£)  ~  -f  ""  U?it  kre"*‘r£(r)] 
*u  2jrircos0  |  r  |  -w  j 

*10$)  =  nj,2  .  AE  ^$)- 

2ak  coeS 

Furthermore, 

*01$)  "  ~T~ 0x9  *M^) 

•o 

—  —5 - i - Ag  i(-k). 

’'o  2*^0580  & 

Substituting,  we  get 


1  Jo. 

4?~  % 


-2*i  aE  1 


abk'  1- 


ttzFoJ 


E„ab  i-  Ur 


,  ..  Y.,cos$  . 

voi<k> 


— -1554g»g - yi-  ^  1(4) 

jbirnxF  3^#*° 


for  an  OEWG,  or 


=S2£ - - - limit  kra~^'r  Efr) 

01  Ajrr^jK^i.1,1  J 

where  E($,<j>)  is  the  transmitted  field  due  to  an  excitation  Efl  in  an  OEWG 
Sex(x)  =  E0cos-^-ey. 

E.  FAR  HELP  QF  PYRAMIDAL  8QRN.MlgM/S,  . 

The  pyramidal  hom  probe  antenna  has  the  advantage  that  it  has  higher  gain,  which 
increases  the  signal-to-noise  ratio  (SNR)  on-axis.  Of  course,  the  horn's  higher  directivity 
reduces  the  SNR  off-axis,  and  the  presence  of  nulls  in  the  hom  pattern  is  undesirable. 
Since  the  hom  antenna  is  well-characterized,  it  is  also  a  good  choice  for  a  test  antenna  for 


near-field  antenna  measurements. 

The  following  discussion  of  the  pattern  of  the  pyramidal  hom  is  taken  from  Balani3 
[6].  The  pyramidal  hom  geometry  is  shown  in  Figure  R-5  (Bakuns'  Figure  12.18).  The 
following  parameters  are  assumed  to  be  known: 

a  s  waveguide  feed  dimension  (larger),  along  ex, 
b  s  waveguide  feed  dimension  (smaller),  along  e  , 
a^  =  hom  dimension  (larger),  along  es, 
bj  s  hom  dimension  (smaller),  along  ey, 

Pji  s  distance  from  aperture  to  E-plane  phase  center, 

P2  5  distance  from  aperture  to  E-plane  phase  cent®, 

Eq  =  waveguide  excitation  (TEjq  mode),  and 
A  5  free-space  wavelength. 

The  tangential  components  of  the  E~  and  H-fields  at  the  hom  aperture  are 
approximately 

.  +ikf-5^  +  -il2] 

Sy(x'df')  =  Eca»[^-*']e  2  1-*. 

_  ,  ,  ~  Eyfc'-T) 

- L—. 


The  far-field  E-field  can  be  derived  (after  tedious  manipulation) 
Ef  =  0, 

E$  =  So  ~Tn“  e+53ff  *inW)  ‘  t1+co*(  *)]  •  Ji*2  1  >  “d 
=  Eo  '4tt'  6+ikr  c06(^),(1+co*(f))‘IiI2]i 

where 


k'=  k  sin (ff)  cos(tf)  +  t/s^. 
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*2- 

k£=  k  sw(S)<xa(<j>)  -  r/ij, 
ky  =  k  sin(0)sin($; 
and 


C(x)  *  cot 

•>0 

L  j 

S(x)  =  sin 

[-#-]* 

•>0 

are  the  Fresnel  cosine  and  sine  integrals. 

We  use  the  Narda  640  horn  antenna  at  10  Ghz  driven  by  X-band  waveguide. 

Typical  parameters  are 

a  =  . 02286  m 
b  =  .01016  m 
aj=  ,0737  m, 

bj=  .0544  m, 

Pj  =  .0918  m, 

P2  =  .1189  m. 

A  =  .03  m. 


APPENDIX  C 

SCATTERING  BY  CONDUCTING  SPHERE 


The  conducting  sphere  is  a  canonical  object  for  which  an  analytical  solution  exists. 
In  this  appendix  we  present  a  solution  [19]  for  scattering  by  a  perfectly  conducting  sphere, 
followed  by  a  derivation  of  the  coordinate  transformations  necessary  to  utilize  this  solution 
in  verifying  the  measured  scattering  matrix. 

For  convenience,  we  will  work  initially  in  an  unprimed  coordinate  system.  A 
perfectly  conducting  sphere  of  radius  a  is  placed  at  the  origin  (Figure  C— 1)  and  is 
illuminated  by  a  monochromatic  incident  plane  wave  traveling  in  the  +e2  direction  with 
x-polarization 


k  =  k  e 
z  z 

.  i(k  JrUt) 

E.  =  E  e  2  e  . 
inc  o  x 

Dropping  the  dependency,  write  £jnj,  in  spherical  form 


where 


-m  yk2R)'pn(C0S^)  sphcosff) 

4/1  \  R)  1 

neln  25  u(a+l)  EIR  ?n(cosff)co8$  ef  + 


S{k2Rjn(k2R)l  aP^coeQ 


cos  lie .  * 


and 


kj  =  propagation  constant  inside  sphere, 
kj  =  propagation  constant  outside  sphere,  and 
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R  =  radius  to  observation  point. 

For  R  >  a,  the  scattered  field  will  have  tb^  form 

n=l 

“eln  "  P(n+t)— e,  + 

1  flP^coafl 

nr - - - - cosp  e.  * 

K2R  ^(kjjR)  88  * 

■ ■*&***  V 


We  define 


fcl-N.k* 

P  =  V** 


kj-a  =  N'/>. 

The  boundary  condition  at  R  =  a  requires  that 
ar  K^l  ~ 

and 

br  ._  ~ 

The  above  equations  are  valid  for  a  sphere  with  arbitrary  dielectric  properties.  If 
the  sphere  is  aluminum,  and  the  surrounding  medium  is  air,  then 


=s  =  1.257  X  lCf*  henry/meter, 

£l  =  <2  =  to  =  8.854  X  10-12  farad/meter, 

<Tj  =  3.5  X 107  mho/meter,  and 
<r2  =  °. 

The  general  form  of  the  propagation  constant  is 
k2  =  i\u£  +  iu/ur 
=  epu^(l  +  i»7)t 

where 

"*-nr- 

At  f  =  10  Ghz, 

’l“s^-S6.3X107 

so 

k2  2  i^/^J2^,  and 

kj  ~  *  JT  •  1.66  X  106  /meter. 

Also. 

•  k2  -  c^u/2 
kj  2  209.4  /meter. 

Choosing  a  6-inch  diameter  sphere,  a  »  .0762  meter,  so 

/>  =  V 

=  15.96 


and 

Np  ~  ±  JT"  •  1.26  X  105  /meter. 

Since  |Np|»l,  we  can  use  the  asymptotic  relations 


and 
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1  1  »-*i»>-°H 

to  simplify  the  above  equations  for  a*  and  b£,  yielding 


r 


b^(/>)sjn(p)  +  i-i?n(p). 

After  the  usual  sort  of  manipulation, 


r,  .  ;  2n+i 

E0  “  Eo~1KJS  c0|v W  . smfl  + 


APPENDIX  D 

SCATTERING  BY  CONDUCTING  DISK 


To  calculate  scattering  from  a  conducting  disk  we  use  a  physical  optics  solution  due 
to  Capt.  Keith  Tiott  of  RADC.  Using  the  theory  of  physical  optics  [20],  on  a  scattering 
object 

J8  a  2  n  *  H 

where 

Js  s  total  surface  current, 

h  s  unit  normal  vector  pointing  away  from  surface, 

§  s  incident  H  field. 


This  is  a  good  approximation  for  large  objects  but  tends  to  break  down  with  smaller 
objects  because  it  is  not  valid  near  edges.  If  the  incident  plane  wave  £  is  known,  then  one 
can  compute,  in  sequence,  H,  3,  and  the  vector  potential 

ff  +£[?-*' | 

A  -  j— -  -i-S - dS' 

**  JJ  i~-r'| 

and  the  scattered  E-field  is 

A. 


In  the  unprimed  frame  the  disk  is  at  the  origin  (Figure  C-l),  lies  in  the  x-y  plane, 


and  has  radius  a.  It  is  illuminated  from  above  the  x-y  plane  by  the  incident  wave  vector 
kj  =  — (cos^sin^-^  +  sin^sin^-ey  +  cosffj’ez)  |k|. 

Assume  kj  lies  in  the  x-s  plane,  so  ^=0,  and 
5.  =  -  (sin A, •  ex  +  0 •  ey  +  cos ^ •  esl  ’k,. 

The  observation  location  is  also  above  the  x-y  plane  at 


r  *  xe  +  y-e  +  a-e 
*  y  * 

=  r  (cos#3in#  e„  +  sin^6in0  e„  +  cos0  e_), 
x  y  z 


Figure  D-l.  Coordinate  System  for  Disk 
Scattering  Prediction 
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k-*r  =  -k.r(8in^cos^6in^  +  cos0gcos$). 

The  plane  of  incidence  is  defined  by  kg  and  eg.  There  are  two  polarizations  of  the  incident 
field  —  parallel  and  perpendicular  to  the  plane  of  incidence.  An  arbitrary  planar 
illumination  EQ  with  kg  (^=0)  can  be  expressed  in  terms  of  those  two  polarizations: 


,  +ik.T 

£||  =  E0j|e  (-coe^  +  sin^) 

8  -E 

±  o±  y 

Considering  the  parallel  polarization  first,  we  have 


H  = 


+ik.-r 
e  e„ 


v  *  ~r 

On  the  disk  ( 8'=x/2 )  we  use  primed  coordinates  to  indicate  a  source,  so 


HI 


1 9=r/2 


E0j|  -ikr'sintfgCos^ 


Then 


jV) 


:2ezxl 


,8=**/  2 
-2E0j[  ’-ikr'sintfgcos^' _ 


The  observation  point  is 

r  =  j  r  |  (cos^sin$ex  +  sindsintfey  +  cosS^) 
while  on  the  disk  itself  0'=jr/2,  so 
?'  =  |r|(cos^'ex  +  sin<S'ey). 

The  factor  |  r-r'  |  in  the  integral  for  A  is 


|r-r'|  = 


■J 


O - TT,  . . .  . . . . . 

(rcos^sintf-r'cos^')  +  (rsmtfsin0-f'sin$')  +  r  cos  0 
r^  +  r'2  -  2rr'8infcos(^-^') 


In  the  denominator  of  the  Integral,  |r— r'|  can  be  approximated  by  |r|.  In  the 
numerator  phase  factor,  a  better  approximation  is  required.  Ignoring  quadratic  terms 
inside  the  radical,  expanding  in  a  Taylor  series,  and  retaining  the  first  two  terms 
| r-r'  |  2  r  -  rsinffcos^- ^'). 
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Then 


.  1 
'  4r?  ' 


f  — 2E, 


'°L 


-4kr 


■e„e 


'shftcos^'  +ik(r-t'8in&oe(^-^/) 

-  e  r'dd' 


d  <p'te' 


After  tedious  manipulation, 

*-«£0|j  -4-ikr  J^(kAi) 
ET 


A. 


Jt 


e„ 


-aE0j|  -S-ikr  J^kAa) 

_i_e 


M — [oafeos^-sin^J 


where 


A  =  |  sin?8  +  sin^J.  +  Zsinitsic&osd  - 
In  the  far  field,  retaining  r“*  terns  only, 

E0=  IjH^a+ikljA^ 

vii)  =  -vn9=+&>i\ 

30  the  scattered  field  due  to  parallel  illumination  Is 

SS||  =  ^olf^nr-1!^)  [-cosfeos^  +  J 

Considering  perpendicular  polarization, 


E 

H  =  -2±- 


Ojl 


cos  5-  e  -sinA  e, 

IX  2  Z 


+ik.-r 

A  * 


+ikr'sinS;COs^ 


Then 


„  2E 
3  = _ 2i 


A  =  ■ 


+ikr'sin&co3& 
n  ey 

E<u  „  Jl(kaA)  . 

— cs^—'kA;  ey 


E  J  (kaA)  ,  , 

-  cosffj  a - gj —  jcosfcin^  +  cos^e^  j 

and  the  scattered  field  due  to  perpendicular  illumination  is 

S8i  =  +iE0Ae+ikr  pjp  Jj(akA)  [cos^cosftin^e^  +  cos^e^)} 

Since  §sj|  and  E^  have  identical  time  phase,  the  scattered  E-fidd  is 


V 

-ikr  , 

=  +i-VirJ  iO*a) 

-cos icoitfi  cos^cosfoin^' 
sin^  cos  A  cos  d1 

m 
.  01 . 
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Mg8§i8»S 
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MISSION 


Rome  Air  Development  Center 


&  RADC  plans  and  executes  research,  development,  test  and  ^ 

?  selected  acquisition  programs  in  support  of  Command,  Control,  ^ 

5  Communications  and  Intelligence  (OI)  activities.  Technical  and  v 

ft  engineering  support  within  areas  of  competence  is  provided  to  k 

*£  ESD  Program  Offices  (POs)  and  other  ESD  elements  to  Q 

^  perform  effective  acquisition  of  (PI  systems.  The  areas  of  f, 

J  technical  competence  include  communications,  command  and  Q 

jg  control,  battle  management  information  processing,  surveillance 

6  sensors,  intelligence  data  collection  and  handling,  solid  state  k 

H  sciences,  electromagnetics,  and  propagation,  and  electronic  q 

f  reliability/ maintainability  and  compatibility.  $ 
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